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ABSTRACT 


A three dimensional finite element based numerieal model was used to analyze 
the heat transfer characteristics of various staggered short pin-fin array heat exchangers. 
The simulation was validated against data from an experimental rig as well as historical 
data and then used to estimate the heat transfer coefficient and pressure drop for a wide 
range of Reynolds numbers for circular and airfoil-shaped pin fins. Circular pin 
configuration variations included changes in pin spacing, axial pitch and pin height ratio. 
Airfoil pin variations also included changes in length and aspect ratio. Correlations for 
Nusselt number and friction factor were developed. Using established performance 
metrics, optimum configurations for both pin shapes were determined. The optimum 
airfoil pin array was shown to match the heat transfer rates obtained by the optimum 
circular pin configuration while incurring less than one third the specific fluid friction 
power loss. The results from this study would be of direct value in the design of a shroud 
enclosed heat exchanger concept being proposed for turbine blade cooling, or for cooling 
of high power electronic components, or in other high heat flux dissipation applications 
requiring a low-profile, high area-density based micro-heat exchanger design. 
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EXECUTIVE SUMMARY 


Motivated by the growing trend towards the development of miniaturized heat 
exehangers for large heat flux applieations, a eomprehensive study was earried out of the 
flow and heat transfer eharaeteristies of staggered pin-fin array heat exehangers of 
various eonfigurations. The primary thrust of this work was based on a fully three 
dimensional finite element based numerieal model of the heat exehangers. The eritieal 
length seale of the heat exehanger that allows sealed up versions of the problem to be 
studied was identified and verified. For a limited set of eonfigurations, the numerieal data 
were validated against experimental data obtained from a modular rig that was designed 
to eover various seenarios. These data were used to prediet the heat transfer eoeffieient 
and pressure drop for a wide range of Reynolds numbers for eireular and airfoil shaped 
fins. Cireular pin eonfiguration variations ineluded ehanges in pin spaeing, axial piteh 
and pin height ratio. Airfoil pin variations also ineluded ehanges in length and aspeet 
ratio. Working eorrelations for Nusselt number and frietion faetor were developed. 
Using established performanee metries, optimum eonfigurations for both pin shapes were 
determined. The optimum airfoil pin array was shown to mateh the heat transfer rates 
obtained by the optimum eireular pin eonfiguration while ineurring less than one third the 
speeifie fluid frietion power loss. The results from this study would be of direet value in 
the design of a shroud enelosed heat exehanger eoneept being proposed for turbine blade 
eooling, or for eooling of high power eleetronie eomponents, or in other high heat flux 
dissipation applieations requiring a low-profde, high area-density based miero-heat 
exehanger design. 
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I. INTRODUCTION 


A, BACKGROUND/MOTIVATION 

Historically, short pin-fin arrays have played an important role in enhaneing heat 
dissipation, partieularly for turbine blade eooling. Their ability to dissipate heat 
effeetively has enabled blade designers to meet the rigorous demands for higher turbine 
inlet temperatures. The requirement for inereased speeifie thrust from future engine 
designs will serve to demand even more performanee from any turbine blade eooling 
scheme. 

Commonly, pin-fin arrays are made up of banks of solid pins that span the end 
walls of an internal flow passage or duet. The array is usually eonfigured with eylindrieal 
pins in either a staggered or in-line arrangement with the eoolant flowing perpendieular 
to the pin axes. Typieal array dimensions are pin diameter (O), pin height {H), 
streamwise pin spaeing (X) and spanwise pin spaeing (S) as shown in Figure 1. Pins ean 
either be long with H/D greater than about 10, or short with H/D on the order of unity. 
Turbine blade eooling arrays fall into the latter eategory and form the subjeet of this 
study. 



Figure 1 Sehematie of a staggered pin-fin array 

Currently turbine blade designers employ these pin-fin arrays mainly in the 
trailing edge of the blade as shown in Figure 2. Due to manufaeturing eonstraints, ribbed 
ehannels and impingement eooling eannot be aeeommodated in these narrow regions 
(Ref 1). One proposed eoneept for enhaneed turbine blade eooling involves mounting 
mieroseale pin fin arrays on the blade surfaee and eovering them with a thin metal skin. 
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Cooling air would circulate from within the blade, through these shrouded arrays and exit 
at various points for film-eooling as shown in Figure 3. 



Figure 2 Typieal application in a turbine blade 


Shroud 



Figure 3 Proposed blade shroud eonfiguration 


2 







Alternatively, these pin fin arrays could be used to cool integrated circuit chips in 
electronic equipment. These components could be mounted directly on top of a 
microscale pin-fin array cooled by forced convection. An economical electro-deposition 
process for manufacturing microscale pin fin heat exchanger arrays is described in by 
Stephens et al. (Ref. 2) 

B. PREVIOUS WORK 

Because of the critical function pin-fin arrays serve, empirical studies dating back 
to the 1980s have been carried out to determine and improve the effectiveness of such 
heat exchangers. VanFossen (Ref. 3) examined staggered pin-fin array configurations for 
H/D values of 0.5 and 2.0. He found that the presence of the pins significantly enhanced 
the array averaged heat transfer coefficient. He also noted that existing long pin 
correlations overestimated the short pin heat transfer coefficients, and thus established the 
need for further research. Sparrow et al. (Ref. 4) later confirmed this finding by studying 
the heat transfer properties for cylinders adjacent to endwalls. They discovered that wall- 
cylinder interactions were confined to one pin diameter from the wall and decreased heat 
transfer compared to regions of the cylinder away from the wall. 

Taking the study further, Metzger, et al. (Ref. 5), examined the streamwise row- 
by-row development of the heat transfer coefficient for a staggered array. Using a 10- 
row array with H/D = 1.0, S/D = 2.5 and X/D = 1.5 and 2.5, they found that the heat 
transfer coefficient peaked between the third and fifth row of the array. Their work 
established a useful baseline empirical model used to validate future experimental results 
such as in the work of Chyu (Ref 6), Chyu & Goldstein (Ref 7), and Li et al. (Ref 8). 

Jubran, et al. (Ref 9) also experimentally investigated the optimal pin-fin spacing 
that would provide the maximum heat transfer rate per unit area. Their results showed the 
optimum X/D and S/D values to be 2.5 for all flow rates tested. However, they used a 
fixed size endwall that meant various pin spacing configurations resulted in different 
numbers of pins. This aspect of their work presents a difficulty in comparing their results 
with other studies such as the work of Metzger, et al. (Ref. 5). In addition, H/D was 
approximately 9.5, and thus more representative of a long pin array. 
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Other experimental studies have focused on alternative pin shapes in an effort to 
improve heat exchanger performance. In 1983, Metzger et al. (Ref. 10) studied heat 
transfer and pressure drop for oblong pin-fins oriented with the major axis parallel to the 
flow direction. They found that both pressure drop and heat transfer rates were lower 
than they were for circular pin-fin arrays. Following this, Arora et al. (Ref 11), studied 
similar oblong configurations with the major axis aligned at various angles of attack 
relative to the flow direction. Arora found that pressure drop was lower and heat transfer 
rates were as good or better for the oblong pins than circular pins when the major axis 
was aligned with the flow. For other angles of attack, the oblong pins offered no 
advantage over the circular pins. Later, Chen et al. (Ref 12) and Li et al. (Ref 8) studied 
drop-shaped and elliptical pin fin arrays respectively. They concluded that these alternate 
shapes provided better heat transfer rates and approximately half of the friction factor of 
the circular pin-fin arrays. 

The relative contribution to overall heat transfer from the pins and the endwalls 
has been an issue of considerable debate, and has been discussed in several earlier 
studies. VanFossen (Ref. 3) first reported that the pin surface heat transfer coefficient 
was 35 percent greater than that of the endwall. Metzger, et al. (Ref 10) determined that 
both were within 10 percent of each other, and then later reported that the pin surfaces 
had approximately 50% higher heat transfer coefficient than the endwall. Contradicting 
all previous results, Al Dabagh, et al. (Ref. 13) indicated that the endwall heat transfer 
coefficient was 15 to 35 percent higher than that for the pin fin s . Offering a solution, 
Chyu, et al. (Ref. 14) determined that the heat transfer coefficient on the pin surface was 
10 to 20 percent higher than the uncovered endwall. 

Because experimental work can be expensive and time consuming, not all 
configurations of interest have been tested. With the current availability of tremendous 
computing resources and several powerful commercially available computational fluid 
dynamics software packages, the next logical step was to conduct these studies 
numerically. Recognizing the potential advantages afforded by numerical modeling for 
pre-design evaluation, Shah et al. (Ref. 15) have made a compelling case for the 
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imperative need for eomprehensive eomputer-based studies of eompaet heat exchanger 
behavior and performance prior to empirical characterization and manufacture. 

Donahoo, et al. (Ref. 16) approached this study as a numerical optimization 
problem using a general-purpose viscous flow solver to simulate the heat transfer and 
flow features within a staggered pin-fin array. A 2-D model was constructed and meshed 
to capture the tumultuous flow characteristics around the pins and in the wall regions of 
the internal flow passages of the heat exchanger. This study has provided some insight 
into the flow and heat transfer characteristics, however due to, the 2-D nature of the 
model it cannot be used to treat the effects of the bounding endwalls and hence the role of 
flow passage height cannot be incorporated into the characterization. 

C. OBJECTIVES 

After an extensive literature search, it was evident that the staggered short pin-fin 
heat exchanger had not been simulated with a three dimensional numerical model. The 
major goal of this research was to develop an experimentally validated three dimensional 
numerical model of the staggered short pin-fin heat exchanger. This model would then be 
used to accomplish the following objectives: 

1. To quantify the heat transfer characteristics and pressure drop of several 
staggered short pin-fin array heat exchanger configurations. 

2. To develop simple and reliable pin-fin heat exchanger performance 
correlations for various circular pin configurations and sizes. Demonstrating that these 
correlations are independent of heat exchanger size would be critical for extending them 
to microscale heat exchangers. 

3. To perform a parametric study of an airfoil-shaped pin fin heat exchanger. 

4. To determine optimal pin-fin array configurations based on suitably 
defined measures of heat exchanger effectiveness. 
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D, METHODOLOGY 


The basic steps involved constructing a suitable numerical model, validating the 
model with experimental results and finally applying the model to a carefully defined test 
matrix. The following sections detail this process. 
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II. NUMERICAL ANALYSIS 

In this section, the numerical portion of the heat exchanger study is presented. 

A, CIRCULAR PIN BASELINE STUDY 

This first subsection provides a detailed numerical analysis of the circular pin-fin 
array with H/D = 1. 

1. Numerical Simulation 
a. Modeling 

The numerical modeling was performed using the FLOTRAN solver of 

the commercial finite element analysis package ANSYS (version 6.0). This code was 

selected with the expectation that periodic unsteadiness would be confined by relatively 

tight physical bounds and would not significantly affect the total time averaged flow 

solution. Based on later observation of the steady reduction in numerical residual values 

and steadiness in experimental data, this assumption was considered to be justified. This 

study examined turbulent flow and heat transfer characteristics within a planar 3-D 

staggered short pin-fin array heat exchanger. 

In order to reduce computational requirements, some simplifying 

assumptions were made. First, by taking advantage of the symmetry planes, only one 

quarter of the heat exchanger was modeled as shown in Figure 4. Further, only the fluid 

(air) flow and heat transfer behavior were modeled thereby eliminating the need to model 

the solid regions of the heat exchanger (pins and end wall) and calculate their temperature 

distributions. This approach assumes that the heat exchanger surfaces could be treated as 

isothermal. Given that the pins are relatively short and metallic, a simple pin-fin analysis 

was carried out to validate this assumption. Calculations justifying this decision are 

shown in Appendix A. Therefore a conjugate analysis has not been implemented at this 

stage, and the passage walls and pin-fins were modeled as solid isothermal boundaries. 

Two basic models were created in this study. The first model, shown in 

Figure 4 simulated a staggered array consisting of 10 rows of pins in the streamwise 

direction, and alternating between 4 and 5 pins per row in the spanwise direction. This 

represented a heat exchanger of finite width similar to the arrays tested by Metzger et al. 

(Ref 5) and was used primarily for validation runs. The second model also simulated a 

staggered pin array with 10 rows of pins, but infinitely wide in spanwise extent. In this 
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Ux = Uy = Uz = 0 Passage Sidewall 

Adiabatic. Ux = Uy = Uz = 0 


Figure 4 Numerical model details 


Flow, T = 300K 



Figure 5 Infinite width numerical model 


8 






















case, one half of an infinitely wide array was simulated. This infinite span eondition was 
numerieally achieved by plaeing “symmetry”, i.e. zero flux boundaries on both the 
spanwise edges of the array as shown in Figure 5. This was done to simulate the 
behavior in mieroseale heat exehangers of ultimate interest, whieh typieally have large 
feature densities eonsisting of many tens of pins per row whereby eontributions from the 
sidewalls would make up but a very small fraetion of the total flow solution. Both the 
finite and infinite models were used to eomplete the ehosen test matrix. Figure 6 shows 
results for a representative eonfiguration that was simulated with both models. It is elear 
that both models produced nearly identical results for Nusselt number and close results 
for friction factor and the implieations will be diseussed in more detail below. Frietion 
factor for arrays with more than 5 pins in the spanwise direetion would fall between the 
limits shown in this graph and are also well within the aeeuraey limits for the frietion 
faetor eorrelation developed in this study. 

For eonsisteney and effieiency, eaeh numerieal simulation was performed 
by exeeuting a maero. Desired test eonfigurations sueh as Reynolds number, axial piteh, 
pins spaeing, height ratio, eooling air and wall temperature were speeified as inputs to the 
maero. In addition to building the model, the maero was used to control mesh density and 
boundary eonditions. Further details on the maero and a sample maero are provided in 
Appendix B. 

To be eonsistent with earlier studies the heat exehanger array test seetion 
was placed between eonstant area adiabatic entrance and exit duct sections (see Figure 5). 
As in Metzger’s experiments (Ref. 5), the entranee duet measured 12.7 em in streamwise 
length. However, while Metzger used a 7.62 em exit duct, the model exit duet in this 
study was extended to 12.5 em, or 25 hydraulie diameters, to ensure well mixed 
eonditions at the exit plane for reliable extraetion of numerieal data. All heated endwall 
and pin surfaees were maintained at a specified temperature. Inlet air veloeity and 
temperature were speeified to produee the desired Reynolds numbers, and the no-slip 
eondition was observed on all rigid boundaries. Veloeities and heat fluxes normal to 
symmetry boundaries were set to zero. 
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3 


X/D = 5, S/D = 3.0 



Reynolds Number 


X/D = 5, S/D = 3.0 



Reynolds Number 


Figure 6 Comparison of results from finite and infinite models 





b. Solution Technique and Numerical Accuracy 

All models were solved using a standardized set of eriteria. Eaoli model 
was solved using the standard kt-z turbulenee model with Van Driest eoupling for the wall 
region. Additionally, all flows were treated as being incompressible and at steady state. 
The CFD FLOTRAN solver in ANSYS was set to use the Preconditioned Generalized 
Minimum Residual (PGMR) method and the Collocated-Galerkin (COFG) approach was 
used to discretize the advection term. More information on these CFD algorithms can be 
found in Appendix C or in the ANSYS Theory Manual (Ref 17). 

Initially, attempts were made to discretize the advection term using the 
Streamline Upwind/Petro-Galerkin (SUPG) approach in order to reduce solution time. 
However, solution analysis revealed that this method although quicker, produced an 
unreliable and oscillatory heat transfer rate solution. Figure 7 depicts the behavior of 
these two solution approaches applied to a typical problem. In this particular case, the 
goal of the task was to replicate the experimental study by Metzger et al. (Ref. 5) in order 
to validate the numerical model. The COFG method yielded a heat transfer rate of 3.41 
watts, while 3.42 watts would have been required to match the results of Metzger. The 
SUPG method appeared to provide a less stable answer, while the COFG method yielded 
stable results that matched experimental results within 0.3%. 

c. Grid Independence 

Meshing was carefully specified to ensure maximum coverage in areas of 
interest such as regions of increased velocity and temperature gradient near the walls as 
shown in Figure 8. Additionally, mesh density was meticulously scaled to facilitate grid 
independence studies. Grid independence was taken to be achieved when the heat 
transfer results could be treated as effectively invariant, i.e. when changes in the overall 
heat transfer rate remained below 2% for subsequent change/refinement in the grid. 
Table 1 provides a representative example of a grid independence check for a given 
configuration. Finally a solution was considered valid only if the heat exchanger fluid 
outlet temperature prediction from the numerical model, matched the temperature 
prediction made from a simple energy balance to within O.IK (<2%). 
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Parameter 

Grid 1 

Grid 2 

Grid 3 

xxyxz 

200x20x9 

220x22x12 

230x23x15 

/pin“ 

m 

/29 

/31 

nodes 

83,760 

123,851 

170,928 

0(W) 

1.8163 

1.9065 

1.9076 

Tout (K) 

303.87 

304.02 

304.01 

m (kg/s) 

4.70E-04 

4.72E-04 

4.73E-04 


“Indicates grid divisions along test section in x, y and z directions 
and along pin edges 
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Figure 8 Sample model meshing 


2. Test Approach 

a. Model Validation 

In order to validate the model, several runs were made in an attempt to 
reproduee the experimental results obtained by Metzger et al. (Ref. 5) and Arora (Ref. 
11). Table 2 shows the eonfigurations and eonditions used for validation. The final 
eolumn shows the temperature differenee between the preseribed isothermal wall and pin 
heat exehanger surfaees, and the ineoming air at 300K. 

For validation, the Reynolds number was defined based on pin diameter, 
and the maximum geometry-based veloeity as follows (originally used by Metzger): 

Re^ = 
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Table 2 Validation configurations 


Config 

S/D 

X/D 

H/D 

Re 

AT(K) 

Metzger 1 

2.5 

1.5 

1.0 

3980, 7310, 13800, 

20000, 30000, 41500, 

69100, 98900 

6,12 

Metzger 2 

2.5 

2.5 

1.0 

3590, 7340, 13900, 

24900 37900, 62100, 

93500 

6,12 

Arora 1 

2.22 

2.83 

1.07 

6500,13827,24000 

6 

Arora 2 

2.42 

2.83 

1.07 

4200,13092,23000 

6 

Arora 3 

2.46 

3.39 

1.28 

4800,11012,29000 

6 


The graphs in Figure 9 through Figure 12 show the results from these 
validation runs plotted and compared with the data and correlations derived from the 
experiments of Metzger and Arora. Here, the Nusselt number and friction factor are 
defined as: 

hD 

= —r 
k 


and 


/ = 


AP 

2pU^ N 

max 


(4) 


where V is the number of pin-fin rows in the streamwise direction. In the 
experiments conducted by Arora, A was replaced with array length. 
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Figure 9 Comparison of numerically obtained Nusselt number with experimental data of 

Metzger (After Ref. 5) 



Figure 10 Comparison of numerically obtained friction factor with experimental data of 

Metzger (After Ref. 5) 
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Reynolds number 

Figure 11 Comparison of numerically obtained Nusselt number with experimental data of 

Arora (After Ref. 11) 


n 



Reynolds number 

Figure 12 Comparison of numerically obtained friction factor with experimental data of 

Arora (After Ref. 11) 
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From these eharts it ean be seen that the numerieal model has predieted 
Nusselt number very well. Table 3 shows a eomparison between the eorrelations 
developed by Metzger et al. (Ref. 5) and eorrelations developed using the numerieal 
model during validation. 

Table 3 Validation results 


X/D 

Metzger experimental 

ANSYS Model 

1.5 

NUj^ 

= 0.069 Re^™ 

=0.075 Re^'"‘ 

2.5 

Nuj^ 

= 0.092 Re^™' 

=0.090 Re^™’ 


The eorroboration with Nusselt number results from Arora is also quite 
good espeeially for the higher Reynolds number eases. However, the numerieally 
predieted frietion faetor results do not provide the same level of eonsistency. Compared 
to Metzger’s data, the numerieal model seems to fall in between the two eonfigurations 
tested. When eompared to Arora’s data, the numerieal model agrees extremely well with 
the experimental eonfiguration 2 but tends to overestimate frietion faetor for 
eonfigurations 1 and 3. Despite these differenees noted, the model was eonsidered useful 
as it displayed relative trends in frietion faetor between eonfigurations to the same degree 
as the experiment. 

Based on the comparisons between the numerical results and existing 
experimental data, the model was considered validated for Nusselt number predictions for 
Reynolds numbers between 3,500 and 100,000. The model was considered useful but 
moderately conservative for predicting friction factors over the same range. 

b. Array Characteristic Length 

For validation, pin diameter was chosen as the characteristic length in 
order to corroborate experimental results. However, with further numerical runs it 
became evident that the pin diameter was not truly representative of the characteristic 
length scale in the convoluted array flow passages for all the test configurations being 
considered in this study since it failed to account for array differences due to variations in 
X/D. Thus a more suitable characteristic length based on the ratio of open volume {Vopen) 
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available for fluid flow in the array, to the total fluid wetted area (Aw), was adopted as 
more truly representative of the tortuous flow paths and defined as: 

4V 

(5) 

A 

This was the same eharaeteristie length used by VanFossen (Ref 3). 

c. Test Matrix 

The model was used to study the staggered short pin-fm heat exehanger 
for various values of pin spaeing and Reynolds numbers. Table 4 shows the variables 
and their values ehosen for this study. Exeept for a few size eomparison runs, all test 
seetions were numerieally construeted with an axial piteh of 12.7 mm, while the pin 
diameter was varied to meet spaeing requirements. In all eases, H/D was set to 1.0. 

Table 4 Test eonfigurations, baseline numerieal eireular pin study 


Variable 

Values 

Reoh 

3,000-50,000 

S/D 

1.25, 1.5,2.0, 3.0, 5.0 

X/D 

1.25, 1.5,2.0, 3.0, 5.0 


Using the revised definition of the eharaeteristie length seale in Eq. (5), 
the Reynolds number is now defined as: 

pUD, 




p 


( 6 ) 


where 


U=^ 


pA 


and 


_ V 

^ _ open 


L 
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3, Data Reduction 

Key results were either obtained direetly from the CFD paekage, or oaieuiated 
using fundamental heat transfer relationships or eonservation laws. CFD direet outputs 
ineluded the total heat transfer rate, inlet and outlet bulk temperatures, pressure drop and 
mass flow rate. The following steps were followed to ealeulate the test array Nusselt 
number and frietion faetor. 

a. Nusselt Number Calculations 

1. Total heat transfer rate and mass flow rate were reeorded direetly from the CFD 
results file after grid independenee had been verified for eaeh run. 

2. Outlet bulk temperature was ealeulated from a simple energy balanee as: 

T =T (7) 

m Cp 

3. Asa eheek, this outlet temperature was eompared with the outlet bulk temperature 
provided by the CFD results file. In most oases, the temperatures were identioal. If 
the temperatures differed by more than 0.1 K, then the run was repeated with either a 
revised/finer mesh and/or with an inoreased number of solver iterations. 


4. With output temperature resolved, the log mean temperature differenoe was 
ealeulated over the test array as follows: 


_ [T^aii -Tbulk ,„) - Tbulk ^^,) 

... f T„,u-Tbulk,„ '\ 


In 


( 8 ) 


5. The heat transfer ooeffioient was ealeulated using the following equation: 


h 


array 


Q 


A 


wetted 


AT, 


Im 


(9) 


19 



Finally, the Nusselt number for the array was determined using: 


Nu^ = ( 10 ) 

b. Friction Factor Calculation 

The frietion faetor in the array was caleulated using the following 

definition: 


/' = 


^array^h 

{PU^L 


( 11 ) 


where L was the overall streamwise length of the array, 
c. Specific Fluid Friction Power 


Speeifie fluid frietion power, useful in eomparing heat exehanger 
performanee was calculated using the following definition: 


E = 


mAp 

pA 


( 12 ) 


where A can be either wetted area or heat exchanger array footprint/face area. 


Note that E has the units of W/m"^. 


4. Results And Discussion 

Upon completion of all the test runs, several key performance indicators were 
examined to understand the heat transfer characteristics and trends for each pin-fin 
configuration and Reynolds number. No data were obtained for the specific case of S/D 
= 1.25, when A/Z) was 3.0 and 5.0, since numerical convergence was not achieved for 
these configurations. 

a. Nusselt Number 

The Nusselt number was calculated for each test configuration. The graphs 

in Figure 13 and Figure 14 show Nusselt number plotted against Reynolds number on a 

logarithmic scale for specific S/D and X/D configurations. While plots of the other 
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configurations are not shown in the interest of brevity, all of them displayed similar 
charaeteristies and trends as the ones provided. 

From these figures it was evident that Nusselt number eould be eorrelated 
as a power law funetion sueh as Equation 13. Moreover, from Figure 13, it was also 
apparent that Nusselt number was fairly insensitive to variations in S/D for a given X/D 
and Reynolds number. However, from Figure 14, it was apparent that S/D had a greater 
effeet on Nusselt number. These observations led to the eonelusion that the leading 
eoeffieient, Cn, was a weak funetion of X/D, and a strong funetion of S/D. Similarly the 
exponent, m, was observed to be a weak funetion of S/D. 

S/D = 2.0 



ReDh 


Figure 13 Nusselt number for S/D = 2.0, cireular pins 

After analyzing the eurve fit data for the different eonfigurations, the 
eorrelation for Nusselt number was derived as: 

(13) 

where the values of Cn and m are provided in Table 5. Flsing these values, 
the array Nusselt number ean be predicted to within 5% for 3,000<Re^^ <50,000 for 

1.25<5'/D<5.0 and 1.25<X/D<5.0. 
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S/D 






1.; 

25 

1, 

,5 

2, 

,0 

3, 

,0 

5, 

,0 

X/D 

Cn 

m 

Cn 

m 

Cn 

m 

Cn 

m 

Cn 

m 

1.25 

0.050 

0.832 

0.069 

0.770 

0.085 

0.748 

0.154 

0.680 

0.193 

0.630 

1.5 

0.071 

0.801 

0.090 

0.753 

0.152 

0.683 

0.166 

0.671 

0.165 

0.650 

2.0 

0.097 

0.777 

0.145 

0.707 

0.183 

0.669 

0.162 

0.680 

0.198 

0.642 

3.0 

- 

- 

0.244 

0.659 

0.252 

0.639 

0.259 

0.626 

0.225 

0.622 

5.0 

- 

- 

0.431 

0.605 

0.210 

0.662 

0.278 

0.616 

0.238 

0.612 


However, these eonstants may also be approximated by simple equations 
for ease in programming applieations using the following definitions: 

Q=0.04(5/Z)f (X/Z^r (14) 
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and 


m = 0M5{S/Dy (15) 


Using these approximations, Nusselt number ean be predieted to within 
±20% for 3,000<Re^^ <50,000 for l.25< S/D<5.0 and 1.25<X/Z)<5.0. 

b. Friction Factor 

The frietion faetor results have also been plotted for eaeh eonfiguration in 
a similar manner. The graph in Figure 15 shows frietion faetor plotted against Reynolds 
number on a logarithmie seale iox X/D = 2.0. As in the ease of Nusselt number, plots of 
the other eonfigurations have been omitted for brevity. It ean be readily seen that the 
frietion faetor varied greatly with ehanges in S/D. This behavior ean be reasoned by 
notieing that for small spanwise spaeing the entering flow is foreed to follow a tortuous 
path around the pin fins thus resulting in large pressure drops. However as S/D is 
inereased there is an effeetive opening up of the flow passages and the frietion faetor 
approaehes the open duet value. The frietion faetor also inereased as X/D was inereased, 
primarily due to the inerease in the overall array length. 

As in the ease with Nusselt number, the frietion faetor was related to 
Reynolds number with a power funetion eorrelation. Using the eonstants provided by 
Table 6, the following eorrelation prediets frietion faetor to within 5% for 
3,000<Reo^ <50,000 for 1.25<<5.0 and 1.25<X/D<5.0. 

= (16) 

Table 6 Frietion faetor eorrelation eonstants, eireular pins, H/D = 1 







S/D 







1.25 


1.5 


2.0 


3.0 


5.0 


X/D 

Cf 

n 

Cf 

n 

Cf 

n 

Cf 

n 

Cf 

n 

1.25 

2.525 

-0.018 

1.273 

-0.040 

1.149 

-0.050 

1.339 

-0.069 

1.281 

-0.155 

1.5 

4.418 

-0.030 

1.849 

-0.045 

1.602 

-0.093 

1.175 

-0.071 

1.152 

-0.140 

2.0 

6.243 

-0.025 

2.318 

-0.041 

1.429 

-0.071 

1.192 

-0.088 

1.189 

-0.141 

3.0 

- 

- 

2.845 

-0.043 

1.538 

-0.071 

1.154 

-0.100 

1.192 

-0.160 

5.0 

- 

- 

3.238 

-0.039 

1.512 

-0.073 

1.155 

-0.119 

1.274 

-0.198 
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X/D = 2.0 



ReDh 

Figure 15 Friction factor ioxX/D = 2.0, circular pins 


As with Nusselt number, there may be cases where a simpler mathematical 
model is desired for programming. In order to derive such a model, relationships for the 
behavior of friction factor as a function of X/D and S/D were examined. Interestingly, 
the behavior of the friction factor for configurations with S/D less than 2.0 was very 
different from configurations with S/D greater than 2.0. Physically, this is likely due to 
the tortuous nature of the flow lines when S/D is less than 2.0 for which the geometry 
clearly shows that the uniform flow at the entrance is forced to negotiate a serpentine 
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path around the pins as noted above. However for S/D greater than 2.0, the geometry 
indieates the presenee of a elear and direet straight “line-of-sight” flow path through the 
array, whereby it ean be argued that some of the flow streamlines ean pass straight 
through the array with minimal or no interaetion with the pins, thus resulting in a 
signifieant deerease in the pressure drop. Figure 16 shows two numerieal solutions with 
streamlines plotted to further illustrate the effeets of S/D on the flow. From this graphie, 
it is elear that the wider pin spaeing results in smoother streamlines and a more direet 
path through the array. The narrow pin spaeing eauses nearly all streamlines to negotiate 
a serpentine path around the pins. These differenees in flow geometry result in very 
different pressure drop eharaeteristies. As a result, two separate eorrelations for frietion 
faetor based on S/D were developed to aeeount for this behavior. 


■vy --'wV- 



^ _ r-, __iTy 


S/D = 5.0 



S/D = 1.5 

Figure 16 Effeet of S/D on streamlines 
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The friction factor correlations developed are as follows: 

/' = ri2-7.6236(5/Z))^°''lRe^°-' + 5(l7) 

where 

B = 98.142 In (X/D) + 0.9197 " (18) 

for 1.25 < SjD < 2.0, and 

B = [0.9262 - 0.1327 (5/Z))] (19) 

for 2.0<5/Z)<5.0. 

These correlations predict friction factor to within ±20% for 
3,000 <Re^^ <50,000 for 1.25 <X/Z)< 5.0. 

c. Heat Transfer Coefficient 

One of the unique advantages of a numerical study is the ability to obtain 
detailed local heat transfer coefficient values at any point in the test array. These data 
provide opportunities to evaluate the relative contributions of both the endwall and pin- 
fin surface to the total heat transfer problem and identify other regions of interest. 
Additionally, it is possible to observe the changes in heat transfer coefficient based on 
location on the pins and along the endwall. Ultimately, these observations can shed light 
on ways to improve heat exchanger performance. In this section, the numerical solution 
for the configuration X/D = S/D = 5.0 at Reoh = 20,000 has been used to probe the local 
heat transfer behavior in the array in detail. 

As previously stated, pin-fins have been introduced in the otherwise planar 
duct in order to enhance overall heat transfer performance. Although endwall surface 
area is reduced by the presence of pins, the pins improve heat transfer rate due to the 
addition of their own surface area, and also by increasing flow turbulence levels on the 
endwall thereby giving rise to better transport rates. Figure 17 provides a contour plot 
indicating the local values of the heat transfer coefficient on the endwall of the array as 
the flow moves from left to right. It is evident that the pins are affecting these local 
values. For example, the heat transfer coefficient is reduced in the wake or “shadow” of 
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each pin where the local velocities are low. Further downstream of each pin, the local 
heat transfer coefficient is enhanced as vortices produced by the pins strike the endwall. 
Note also that as expected the large velocities achieved in the flow region between the 
pins tend to increase the film coefficient. 



Legend: Heat Transfer Coefficient (W/m^k) 

100 200 300 400 500 600 700 

Figure 17 Endwall heat transfer coefficient 

In order to quantify the heat transfer benefits gained by adding pins to the 
planar duct, a comparison was made with empty duct performance. To this end, a 
numerical model with the same dimensions, but without any pins, was constructed. 
Using this model, simulations were performed at Reoh = 20,000, the same conditions that 
were used to develop the contour map in Figure 17. As a cross-check, heat transfer 
performance of the empty duct was also calculated using well established semi-empirical 
correlations (Ref. 18). The numerical and analytical solutions for the empty duct were 
found to be in very good agreement as shown in Table 7. 

Table 7 Empty duct solution comparison 


Parameter 

Analytical 

Numerical 


Prediction 

Solution 

Eriction factor 

0.026 

0.027 

Nusselt Number 

50.5 

53.6 


Eigure 18 provides a qualitative comparison of the endwall heat transfer 

coefficients for the pin-fin array and the empty duct. The data for this plot were obtained 

directly from the resolved numerical model by recording heat transfer coefficient values 

for nodes adjacent to the endwall. These nodes were situated on three separate 
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streamwise paths for the pinned array. One path ran unobstrueted between pins while the 
other two ran through pin eenterlines. These paths were ehosen to qualitatively evaluate 
the effeets of the pins on endwall heat transfer eoeffieient. For the empty duet, a single 
streamwise path was ehosen in the eenter of the duet. Note that for both eonfigurations, 
loeal heat transfer eoeffieients deerease with streamwise travel, initially due to the 
thiekening boundary layer and developing effeets, and after the first row additionally due 
to flow separation effeets. However, in nearly all eases, these eoeffieients for the pinned 
endwall are greater than those on the endwall of the empty duet. This elearly illustrates 
the benefieial effeet that turbulenee from the pins has on the endwall heat transfer 
performance. The only regions of lower heat transfer coefficient for the pinned array are 
in the immediate wake or “shadow” of the pins where local velocities are lower. Slightly 
further downstream of the pin, the heat transfer coefficient rapidly climbs almost to the 
value of the unobstructed path. 



Figure 18 Endwall heat transfer coefficient trends 
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Examination of heat transfer coefficient trends on the pin surfaces also 
provided interesting results. Figure 19 consists of radial plots showing the angular 
variation of the heat transfer coefficient on pins from various rows. In each plot, the flow 
is from left to right. Local heat transfer coefficients were recorded circumferentially at 
two heights along the pin measured from the endwall, at z = D/8 and z = D/2. From these 
plots it is clear that local heat transfer coefficient is highest on the leading face of the pin. 
This region is defined on each plot by angles between 90 and 270 degrees. Continuing 
beyond this region, heat transfer coefficient drops off dramatically as the flow separates 
from the pin surface. Additionally, the values at the z = D/8 height location are 
significantly lower than those at z = D/2 due to the boundary layer effects from the 
endwall. 


Heat Transfer Coefficient Profile on Pin in Row 1 


Heat Transfer Coefficient Profile on Pin in Row 2 




Heat Transfer Coefficient Profile on Pin in Row 4 


Heat Transfer Coefficient Profile on Pin in Row 10 




Figure 19 Circular pin heat transfer coefficient profiles 

The effects of the endwall are further illustrated by plotting average heat 
transfer coefficient development as a function of height along the pin. Figure 20 shows 
this development for pins in various rows. As noted before, the heat transfer coefficient 
increases with distance from the endwall and shows a slight increasing trend even at the 
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mid-section of the duct, i.e. at z = D/2. This is consistent with the findings of Sparrow et 
al (Ref 4) who studied a long isolated cylinder in cross flow to investigate edge effects 
due to the bounding planes of the surface at the ends of the cylinder. As mentioned 
previously, they concluded that wall-cylinder interactions were confined to within one 
diameter of the wall and decreased heat transfer on the cylinder near the wall when 
compared to regions of the cylinder away from the wall. 
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Figure 20 Circular pin heat transfer coefficient trends 

Another interesting observation is that the pins in the second row have the 
highest heat transfer coefficients when compared to pins in the other rows of the array. 
There are two likely contributors to this observation. First and primarily, the second row 
of pins encounters a higher local velocity than the pins in row 1. This is due to the 
“nozzle effect” created by the first row of pins as the flow accelerates to provide a 
constant mass flow rate through a reduced cross sectional area. This accelerated flow 
directly impinges on the second row pins due to the staggered geometry thus significantly 
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increasing heat transfer due to stagnation flow characteristics. A weaker secondary 
reason is that the turbulence from the first row enhances the heat transfer of the second 
row. Beyond the second row however there are no additional increases in velocity or 
turbulence level as shown in Figure 21 and Figure 22 since portions of the flow are now 
separated and do not lend to the same stagnation flow type enhancement. 



U*= U/U. 

m 

0 0.5 1 1.5 2 

Figure 21 Velocity profile, circular pins 



0 0.02 0.04 0.06 0.08 0.10 0.12 0.14 

Figure 22 Turbulence level, circular pins 

Also of interest is the comparison between the pin and endwall heat 
transfer coefficients which has been a topic of much discussion in earlier studies. Figure 
23 shows a representative case. Note that in this case the heat transfer coefficients from 
the pins were clearly about 100 - 200% higher than for the endwall. While this 
magnitude varied for each test configuration, the pin heat transfer coefficient values were 
always found to be higher by about 50-200%. 
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Figure 23 Heat transfer eoeffieient, eireular pins vs. wall 

This eonelusion is eonsistent with a majority of the earlier studies that 
have demonstrated a greater heat transfer eontribution from the pin surfaees. 

d. Optimization 

One of the goals of this study was to prediet optimal heat exehanger 
design eonfigurations that would maximize heat transfer rate while minimizing frietional 
losses in the flow. As reeommended by Kays and London, (Ref. 19) various heat 
exehanger eonfigurations ean be readily evaluated and eompared by plotting the heat 
transfer eoeffieient versus frietion power on a suitable unit surfaee area basis. Sueh a plot 
is provided in Figure 24. For elarity, only two eonfigurations showing the extremes of 
heat exchanger performance are shown, and compared with empty duct performance. All 
other configurations fall between these extremes. From this chart, it was clear that for all 
of the configurations tested, the case of X/D = 5, S/D =3 provided the greatest heat 
transfer rate per friction power while the case of X/D = 1.25, S/D = 5 was the least 
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desirable. Surprisingly, the empty duet provided performanee nearly as good as the least 
favorable pin-fin eonfiguration at higher Reynolds numbers. This highlights the need for 
eareful heat exehanger design in order to verify and fully realize the benefits of using pin- 
fins for performanee enhaneement. 



Figure 24 Performance comparison, circular pins 


It appeared that the relative goodness of each configuration was closely 
related to the volumetric density of available heat transfer surface area, or A^/V. Of all 
configurations tested, the optimum array had the third highest area density of 829 m /m 
while the worst configuration had the lowest density of 222 m /m . 

The significance in the difference in performance of the different 
configurations can be quantified by choosing a specific heat transfer coefficient, and 
comparing the friction power required for each configuration. Using this measure it was 
found that the least efficient configuration requires up to 15-20 times more friction power 
for the same heat transfer coefficient, i.e. relatively small gains in heat transfer coefficient 
require large increases in friction power. For example, for the best configuration, 
increasing the heat transfer coefficient from 300 to 450 W/m K results in friction power 
increases from 3 to 22 W/m^. In other words, increasing heat transfer coefficient by 50% 
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requires over 600% more fluid power. Careful examination of the slope of these graphs 
can help the designer choose an operating region that provides the best cost/benefit ratio. 


e. Application 

The actual performance of these heat exchangers can best be quantified by 
inserting sample initial conditions and calculating predicted performance parameters. For 
example, assuming a 30K difference between wall and cooling air temperature, the 
optimum configuration transfers 0.7 W/cm^ at Reoh = 3,000 and 4.0 W/cm^ at Reoh 
=50,000. 

2 

While 4.0 W/cm is fairly respectable, it is expected that higher rates 
would be attainable by further increasing the area density of these heat exchangers. Two 
methods for increasing the heat transfer surface area density are reducing the heat 
exchanger size and reducing H/D. Reducing either size or H/D by a factor of 2 results in 
a doubling of area density. Numerical studies to quantity the effects of variable size and 
H/D on pin-fin heat exchanger performance were completed and are presented in the next 
sections. 

5, Conclusions 

The objectives of this investigation were met by simulating various pin-fin heat 
exchanger configurations with an experimentally validated three-dimensional finite 
element model. Specifically: 

1. Comprehensive power function correlations were developed for Nusselt 
number and friction factor accurate to within 5% for all configurations at Reynolds 
numbers between 3,000 and 50,000. 

2. Heat transfer coefficients were found to be 50-200% greater for the pin 
surfaces than for the endwalls. 

3. In terms of the heat transfer coefficient, the configuration that required the 
least amount of fluid friction power was X/D = 5.0, S/D = 3.0 Conversely, X/D = 1.25, 
S/D = 5.0 had the greatest pressure losses. 

4. In all pin-fin array configurations, it was found that small increases in heat 
transfer coefficient resulted in disproportionately large increases in frictional losses thus 

suggesting the presence of a sensitive design operating point. 
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B, CIRCULAR PINS, EFFECTS OF AXIAL PITCH 

In the previous seetion, the study foeused entirely on pin-fin arrays with a fixed 
axial piteh of 12.7 mm. It was briefly mentioned that redueing axial piteh would inerease 
area density and therefore had the potential to improve heat exehanger performanee. 
Additionally, it is eritieal to understand the effeets or axial piteh reduetion in order to 
extend the results from maeroseale experiments to mieroseale eonfigurations whieh are of 
ultimate interest. 

1. Test Approach 
a. Test Plan 

The previously defined numerieal model was used to study the staggered 
short pin-fin heat exehanger for various configurations and Reynolds numbers. Table 8 
shows the variables and their values chosen for this study. 

Table 8 Axial pitch test matrix, circular pins 


X/D 

S/D 

H/D 

X 

Reoh 

2 

2 

1 

3-50mm 

3,000-20,000 

2 

2 

1 

0.75-6 mm 

50-2,000 

1.5 

1.5 

1 

3-50mm 

3,000-20,000 

1.5 

1.5 

1 

0.75-6 mm 

50-2,000 


b. Laminar Modeling 

As mentioned previously, the majority of runs were conducted using the 
standard k-e turbulence model. However, since some of the variable axial pitch runs 
included low Reynolds number, laminar modeling was also considered. Figure 25 shows 
an h vs. E chart for X/D =1.5, X/D = 1.5 at X = 3mm. This chart includes numerical 
results obtained with turbulent modeling for Reynolds numbers ranging from 50 to 
20,000 as well as numerical results using laminar modeling for Reynolds numbers 
ranging from 50 to 2000. The models produced nearly identical results for Reynolds 
numbers in the 1,000 to 2,000 range. However, below a Reynolds number of 1000, the 
laminar and turbulent solutions were no longer identical. As Reynolds number was 
decreased, the difference between the two solutions was seen to increase. 

According to ANSYS documentation (Ref. 17) the decision on which type 

of modeling to use can be made by considering the ratio of effective viscosity and 
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laminar viscosity. In cases where the effective viscosity is less than 5 times greater than 
laminar viseosity, the doeumentation reeommends using a laminar solution. In this 
applieation, it was seen that this point occurred near Reynolds number equal to 1000. 
This fact, coupled with the behavior observed in Figure 25, lead to the conclusion that 
laminar modeling should be used for Reynolds number below 1000. 



Figure 25 Comparison of laminar and turbulent models 

2. Results and Discussion 

a. Effects of Axial Pitch on Nusselt Number and Friction Factor 
As shown in Figure 26, the dimensionless quantity Nusselt number 
showed no significant change as axial pitch was varied from the baseline value of 12.7 
mm. As axial pitch was decreased, pin diameter had to be deereased to maintain the 
desired X/D values. As a result, the inlet velocity was increased to maintain the Reynolds 
number at these smaller hydraulic diameter values. This resulted in larger heat transfer 
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coefficients. However, this inerease was offset by the smaller hydraulie diameter when 
ealeulating Nusselt numbers. 

The effeet of axial pitch on friction factor is shown in Figure 27. While 
there are perturbations for the lower Reynolds numbers, there are no indieations that 
frietion faetor is dependent on axial pitch. These perturbations are likely aberrations in 
the numerieal data that are consistent with the level of diffieulty achieving eonvergenee 
for these low Reynolds number cases. 

While these trends were expeeted, it was benefieial to eonfirm them 
numerically. The obvious benefit of these results is that data obtained either numerically 
or experimentally for large-scale models may be applied to mieroscale arrays. 
Additionally, the eonsistent behavior of the dimensionless qualities further eonfirmed the 
validity of the numerical model and its ability to provide aeeurate solutions over a wide 
range of absolute flow veloeities. 





Axial Pitch (mm) 


Figure 26 Effeet of axial pitch on Nusselt number, circular pins 
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Axial Pitch (mm) 

Figure 27 Effect of axial pitch on friction factor, circular pins 


b. Performance Comparisons 

As in the previous section, a chart plotting heat transfer coefficient versus 
fluid friction power per wetted area can be made to evaluate heat exchanger performance. 
Here the desired goal is to maximize heat transfer rate while minimizing losses due to 
pressure drop. Figure 28 depicts h vs. E for the configuration X/Z) = 1.5, S/D = 1.5 at 
various axial pitches and Reynolds numbers. The solid lines indicate constant Reynolds 
number with varying axial pitch. Conversely, the dashed lines indicate constant axial 
pitch with varying Reynolds number. From this figure it is clear that heat transfer 
coefficient and fluid friction power increase as axial pitch decreases and Reynolds 
number increases. Further, it can be seen that in both all cases, relatively small increases 
in heat transfer coefficient result in large increases in fluid friction power. A similar 
graph can be prepared for XD = 2.0, S/D = 2.0 but has been omitted for brevity. 


38 











X/D = S/D = 1.5 



Friction Power/Wetted Area (W/rrf) 

Figure 28 Effect of axial pitch on heat exchanger performance, circular pins 

Further examination of Figure 28 reveals that it is more “cost effective” to 
enhance heat transfer coefficient through reductions in axial pitch rather than through 
increases in Reynolds number. This is mainly due to the advantages of scale on both heat 
transfer coefficient and fluid friction power. When holding Reynolds number constant, it 
was observed that heat transfer coefficient varied inversely with axial pitch. However, 
because pin diameter was directly related to axial pitch for a given configuration, fluid 
friction power was inversely proportional to axial pitch cubed. Combining these trends 
produces the following relationship for heat transfer coefficient and fluid friction power. 

h oc (20) 

When holding axial pitch constant, it was observed that heat transfer 
coefficient varied with the square root of Reynolds number. However, fluid friction 
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power varied with the Reynolds number squared. Combining these trends resulted in the 
following relationship. 


h oc (21) 

Both of the above relationships are eonsistent with the trends observed in 
Figure 28. While the differenees in the exponents appear to be small, the impaet of these 
differenees is very signifieant. By way of example, eonsider the line of eonstant 
Reynolds number =100 shown in Figure 28. The furthest point to the right has an axial 
piteh of 0.75 mm and produees a heat transfer eoeffieient of approximately 370 W/m K 
at a frietion power eost of about 0.2 W/m^. To aehieve an equivalent heat transfer 
eoeffieient of this mieroseale heat exehanger using an array with an axial piteh of 3 mm, 
a Reynolds number of about 2,000 is required whieh produees a frietion power loss of 
approximately 6 W/m . Thus, seleetion of an array that is only 4 times larger results in 
an inerease of fluid frietion losses by a faetor of 30. 

3. Conclusions 

The objeetives of this phase of the investigation were met by simulating pin-fin 
heat exehanger eonfigurations with various axial piteh values with an experimentally 
validated three-dimensional finite element model. Speeifieally it was found that: 

1. Variations in axial piteh were shown to have no appreeiable effeet on 
Nusselt number and frietion faetor. This not only validated the choiee of hydraulie 
diameter as a suitable eharaeteristie length but also demonstrated that the results of 
maeroseale experiments eould be direetly applied to mieroseale heat exehangers. 

2. Reduetions in axial piteh ean produee a signifieant inerease in heat 
exehanger performanee. In faet, it was shown that redueing axial piteh is more eost 
effeetive than inereasing Reynolds number to improve performanee. 
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C. CIRCULAR PINS, EFFECTS OF PIN HEIGHT RATIO 

The purpose of this seetion is to discuss the effect of pin height ratio on circular 
pin-fin heat exchanger performance. Variations in pin height ratio, like reductions in 
axial pitch have a large effect on area density. This phase of the study demonstrates the 
result of these effects. 

1. Test Approach 

The majority of test runs were conducted using a numerical model with pin 
spacing ratio of S/D = X/D =2 and H/D varying from 0.25 to 4.0. Reynolds number 
ranged from 2,000 to 64,000. Additional configurations were also run to spot check 
Nusselt number and friction factor correlations developed from this primary test 
configurations and are shown later. 

2. Results and Discussion 

Upon completion of the test runs, key data were obtained and Nusselt number and 
friction factor values were calculated for each configuration and Reynolds number. The 
results of these runs are presented in this section. 

a. Effects of Pin Height Ratio on Nusselt Number 

Figure 29 shows Nusselt number plotted against Reynolds number for 
various H/D configurations of the single pin spacing X/D = S/D = 2.0. Note that the 
Nusselt number tended to increase as H/D increased from 0.25 to 4.0. This effect is more 
pronounced at higher Reynolds number. 

This increase was consistent with the geometric differences in the test 
arrays. As H/D increased, hydraulic diameter also increased resulting in lower absolute 
flow velocities for a given Reynolds number. These lower velocities produced lower 
array averaged heat transfer coefficients. However, the hydraulic diameter increased at a 
greater rate than the heat transfer coefficient dropped. As a result, the Nusselt number 
increased. 
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X/D = S/D = 2.0 





ReDh 


Figure 29 Effect of H/D on Nusselt number, circular pins 


Compared to the baseline case of H/D = 1.0, each variation in H/D acted 
as a magnification factor that depended on Reynolds number and pin height ratio. It was 
recognized that this amplification factor could be estimated by a power function 
correlation. By curve fitting these factors, the following magnification relation was 
derived; 



( 22 ) 


Combining this amplification factor with the previously developed Nusselt 
power function correlation of equations 13 through 15, results in: 





( 23 ) 


where 
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and 


m = 0.865 
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This correlation predicts Nusselt number to within ±30% for 3,000 < < 50,000 

for 1.25<-|<5.0, l. 25 <f < 5 .oand o.25<f <4.0 


b. Effects of Pin Height Ratio on Friction Factor 
Figure 30 shows friction factor plotted against Reynolds number for all 
H/D configurations tested. As with Nusselt number, frietion factor is magnified as H/D 
increases from 0.25 to 4.0. 

As with Nusselt number, this inerease was also consistent with geometric 
characteristics of the test arrays. Higher H/D resulted in lower flow velocities and 
subsequently less pressure drop for a given Reynolds number. However, the rate of 
deerease in pressure drop was overeome by the rate of deerease in the square of the 
velocity. As a result the friction factor increased for higher H/D configurations even 
though pressure loss was dramatically lower than in the lower H/D configurations. 
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Figure 30 H/D effects on friction factor, circular pins 
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As before, it was apparent that a power funetion eorrelation would 
adequately represent the magnifieation of frietion faetor. Curve fitting these faetors 
resulted in the following amplifieation equation: 


A. = exp 


r 


0.2 


H 


In-ln2 

V D 


0.1191— -1| 

Re^ > (24) 


Combining this amplifieation faetor with the previously developed frietion 
faetor power funetion eorrelations results in: 


/' = Re;,^ 

where 


(25) 


Cp = 22.67 
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+ 1.4 


and 


n = -0.085 
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for 1.25 < — <2.0 
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and 


C, = 0.25 
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+ 1.7 


and 


n =-0.041 
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kDj 


-0.019 


for 2.0 < — <5.0 
D 

These eorrelations prediet frietion faetor to within ±30% for 
3,000<Re^^ <50,000 for l.25<f <5.oand o.25<f <4.0. 

Since these correlations were developed based on trends seen in a single 
configuration, it was necessary to validate its ability to predict Nusselt number and 
friction factor over a wide range of configurations and Reynolds numbers. To 
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accomplish this, the conhgurations listed in Table 9 were run. The last two eolumns 
indieate ratios of values predieted by the eorrelations divided by values obtained 
numerically. In the majority of oases, the values predioted by the correlations are within 
30% of the numerically obtained values. 

Table 9 Correlation “Spot check” results, circular pins 


X/D 

S/D 

H/D 

Reoh 

NUc/NUm 

fc/fm 

1.5 

1.5 



1.01 

1.02 

1.5 

1.5 



1.34 

0.67 

1.5 

1.5 



0.95 

1.14 

1.5 




0.95 

0.95 

1.5 




1.33 

0.65 

1.5 




0.79 

0.96 

1.5 




1.00 

0.75 

3.0 




1.11 

0.75 

3.0 




0.92 

0.79 

5.0 

1.5 



0.69 

0.73 

5.0 

1.5 



1.02 

0.87 

5.0 




0.95 

0.99 

5.0 




0.73 

0.85 

5.0 




0.97 

0.58 


c. Combined Effects of Axial Pitch and Pin Height Ratio 

One oritioal objective of this study was to ensure that heat exchanger 
performance correlations were applicable to microsoale designs. Earlier, it was shown 
that Nusselt number and friotion factor were relatively independent of axial pitch at H/D 
= 1. A spot check was run at S/D = X/D = 5.0 and H/D = 4.0 to demonstrate that these 
important dimensionless parameters remained unehanged as axial pitch was varied. 
Figure 31 and Figure 32 show the results of this investigation. It was observed that 
Nusselt number was insensitive to ehanges in axial piteh and friotion faotor ohanged by 
less than 15% as axial pitch was varied from 3 to 50 mm. This Ending is important if the 
previously developed correlations are to be used to predict heat exchanger performance in 
microsoale applioations. 


45 








Axial pitch (mm) 


Figure 31 Effect of axial pitch on Nusselt number, H/D = 4 



Axial pitch (mm) 

Figure 32 Effect of axial pitch on friction factor, H/D = 4 





d. Performance Comparisons 

In order to determine the effeets of variable H/D on heat exehanger 
performanee, all eonfigurations tested were plotted as a series of h vs. E eurves. Figure 
33 eontains a plot of h vs. £■ on a wetted area basis while Figure 34 shows h vs. £■ on a 
face area basis. On a wetted area basis, heat exchanger performance tended to decrease 
as H/D was increased from 0.25 to 4.0. This performance decline was mainly due to the 
decrease in area density as pin height ratio increased. For H/D = 0.25, the area density 
was 11.50 cm /cm , while for H/D = 4.0 it decreased to 1.75 cm /cm . 



Friction Power/Wetted Area (W/mP) 

Figure 33 Effect of H/D on performance, wetted area basis 


However, when considering performance based on a face area basis the 
optimum pin height ratio tested was 4.0 despite the fact that this configuration had the 
lowest area density. This was likely due to two factors. First, the longer pins presented 
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Optimization 



Figure 34 Effect of H/D on performance, face area basis 

increased surface area to provide heat transfer. Thus a greater amount of heat had to be 
transferred through the heat exchanger footprint thereby increasing the effective heat 
transfer coefficient. Additionally, it has been shown Sparrow et al (Ref. 4) that the 
presence of the endwall tends to degrade heat transfer coefficient on the pins. As this 
degrading effect is primarily limited to a distance from the endwall equal to 
approximately one pin diameter, the longer pins would have greater length with a higher 
heat transfer coefficient. This point is further illustrated in Figure 35 which shows the 
heat transfer coefficient for pins of varying height ratio. Note that for H/D = 1, only a 
small percentage of the pin length is close to the maximum value of the mid-pin heat 
transfer coefficient. As pin height ratio increases, a larger percentage of the pin length is 
exposed to maximum heat transfer levels. Note that H/D = 8 was not part of the test 
matrix but it has been included in this chart to further illustrate the trend. 
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X/D = S/D = 5.0, Constant Inlet Velocity 



Figure 35 Heat transfer eoeffieient development on eireular pins 

Returning to Figure 34, some interesting trends are noted regarding pin 
height ratio effeets. As expeeted from the previous argument, the eonfigurations H/D = 
2.0 and 1.0 provide eonsiderably less desirable performanee than H/D =4.0. However, 
when pin height ratio is further redueed from 0.5 to 0.25, the trend aetually reverses. In 
this ease the inerease in area density assoeiated with the smaller pin height ratio has 
overeome the degradation due to wall interaction on the cylinders. 

While these trends were observed for X/D and S/D = 2.0, it was not 

obvious that the same would hold true for other pin spacing such as X/D = S/D = 5.0. 

This configuration was investigated and the results are provided in Figure 36. In this 

case, the performance of the larger pin height ratio is not significantly different from the 

configuration H/D = 1.0. As before, the changes in area density are likely the reason for 

these trends. As shown in Figure 37, this difference in area density is more significant 

49 




for the widely spaeed pin eonfiguration. ForX/Z) = S/D = 5.0, the area density at H/D = 1 
is signifieantly higher than it is at H/D = 4. As a result, the gains in effeetive heat 
transfer performance achieved through longer pins are matched by the performance gains 
associated with the higher area density of the short pin configuration. Based on these 
observations, it appears that enhancing heat transfer rates across the face area of pin-fin 
heat exchangers by increasing H/D is more effective for streamwise and spanwise pin 
spacing ratios less than 5.0. 



Figure 36 H/D effects on performance, wide circular pin spacing 
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Figure 37 H/D and pin spacing effects on area density 

e. Effects of Pin Height Ratio on Endwall Heat Transfer 
Coefficient 

It has been shown above that the presence of the endwall has a significant 
effect on the pin heat transfer coefficient profile. The effect of the pin height ratio on the 
endwall heat transfer coefficient can also be examined using standard CFD post 
processing tools. Figure 38 shows contour plots indicating endwall heat transfer 
coefficient values for three different pin height ratios. Note that H/D = 8 was not part of 
the original test matrix but has been included to explore possible trends. From this figure 
it appears that higher pin height ratios result in greater variations in heat transfer 
coefficient along the endwall especially in the spanwise direction. This is likely due to 
the greater space available for vortex growth as the cooling air swirls around the longer 
pins. Additionally, the impact of viscous damping on the wall may be reduced in the 
longer pin configuration. 
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H/D= 1 



Legend: Heat Transfer Coefficient (W/m^k) 


0 200 400 600 

100 300 500 

Figure 38 H/D effects on endwall heat transfer coefficient 

However, despite the significant local variations in heat transfer 
coefficient seen in long pin configurations, the overall values are not significantly 
different from the short pin case. Figure 39 shows a comparison of average heat 
exchanger coefficient values for different pin height ratios. Cooling air inlet velocity was 
held constant for all cases to provide an opportunity to compare absolute values of heat 
transfer coefficient. A comparison is also made with the empty duct to further illustrate 
the beneficial effect of adding pins to enhance heat transfer rate. Interestingly, the long 
pin heat transfer coefficient profile tends to behave like the profile seen in the empty 
duct. Although its average is significantly higher, the long pin heat transfer coefficient 
profde rapidly decreases and then levels out much like the empty duct case. Also, the 
large variations seen in the contour plots are seen as highly fluctuating absolute values in 
Figure 39. These actions support the earlier suggestion that longer pin configurations 
have reduced effect from viscous wall damping. 

The three dimensional model provides visibility into the nature of the flow 
in the vicinity of the junction of the pins and endwall. Figure 40 shows an isometric and 
side view of streamlines in the flow as they pass over the pins. Note that there is some 
variation of the streamlines in the “z” or pin axis direction. Presumably this adds to the 
effective transport of heat that is promoted by the presence of the pins. Additionally, the 
streamlines are seen spreading around the pin bases in the vicinity of the wall as the 
vortices are formed. 
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f. Optimization 

By combining the effects observed in this study, an optimum theoretical 
heat exchanger configuration can be determined. The optimum configuration would 
consist of small axial pitch, small H/D, large S/D and X/D and operate at low Reynolds 
number. Table 10 shows the logical progression as these enhancing features are applied 
starting from a baseline configuration. All configurations compared have S/D and X/D 
held at 5.0. Within the bounds of this test program, the optimum configuration would be 
defined as X= 0.75 mm, H/D = 0.25, X/D = 5.0 and S/D = 5.0 and would have an area 
density of 525 cm /cm . 

Table 10 Optimization configurations, circular pin-fin array 


Config 

A (mm) 

H/D 

AJV 

(cm^/cm^^ 

baseline 

12.7 

1.0 

8.1 

min pitch 

0.75 

1.0 

137.1 

min H/D 

0.75 

0.25 

525.8 


Figure 41 shows an h vs. E chart that shows the performance of these 
theoretical optimization designs. As expected the baseline configurations provides the 
worst performance. Reducing the axial pitch greatly increases the heat exchanger 
performance. Finally, reducing the pin height ratio results in the best heat exchanger 
performance. Not surprisingly, the performance of these arrays closely follows area 
density. 

The degree of performance enhancement can be appreciated by comparing 
the optimum configuration with the baseline case. For example, at a Reynolds number of 
only 50, the optimum configuration produced a heat transfer coefficient of nearly 1000 
W/m K. The baseline configuration required a Reynolds number of 20,000 to provide the 
same heat transfer performance at a cost of about 300 times the specific fluid friction 
power. For a 3 OK temperature difference between cooling air and heat exchanger 
surface, the optimum configuration would transfer 3 W/cm^ at Re£)/,=50 and 17 W/cm^ at 
R^Dh ~ 1000 . 
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X/D = S/D = 5.0 


d 



Fluid Friction Power/Wetted Area (W/rrf) 

Figure 41 Optimization configuration performance, circular pins 

However, the optimum theoretical configuration would have an actual 
passage height of only 38 pm and would be highly susceptible to particulate 
contamination. Typically, turbine engine designers use minimum passage diameters of 8 
- 12 mils or 0.2 to 0.3 mm to reduce the risk of blockage. Even the intermediate 
configuration, H/D = \,X= 0.75 mm would only have a passage height of 150 pm. With 
this in mind, the theoretical optimum configuration can be modified to achieve required 
clearance. With S/D and X/D held at 5.0, many practical configurations can be defined 
that all have a minimum clearance value of 0.3 mm. Table 11 shows three of these 
configurations and their associated area density. These configurations were chosen to 
take advantage of reducing pin height ratio, axial pitch or a balance of both factors. 
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Table 11 Practical optimum configuration options, circular pin-fin array 


Conf# 

X(mm) 

H/D 

Aj'V (cm"^/cm^^ 

PI 

0.75 

2.0 

73.0 

P2 

3.0 

0.5 

66.7 

P3 

6.0 

0.25 

65.6 


The performance results for these practical configurations are shown in 
Figure 42. From this chart it is clear that P3 was the best performer among the practical 
heat exchanger configurations. 


X/D = S/D = 5.0 



Fluid Power/Aw (W/mP) 

Figure 42 Performance comparison, practical configurations, circular pins 

Figure 43 shows a comparison between the best practical configuration 
and the previously discussed theoretical optimum configuration. The baseline and 
intermediate configurations are included for reference. It is important to recall that 
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neither the optimum nor the intermediate eonfigurations were feasible due to 
suseeptibility to partieulate eontamination. 

Despite the size requirement, this optimum praetieal eonfiguration 
provided exeellent performanee from Reynolds numbers ranging from 50 to 2000. On 
the low end, this eonfiguration had an average heat transfer eoeffieient of 220 W/m K 
and required 0.08 W/m frietion power. At the high end, this eonfiguration produeed an 
average heat transfer eoeffieient of 1050 W/m K and ineurred frietional losses of 43 
W/m . On average the optimum praetieal configuration required over 40 times less 
specific friction power to achieve the same heat transfer coefficient values as the baseline 
configuration. 


X/D = S/D = 5.0 



Fluid Friction Power/Wetted Area (W/mP) 

Figure 43 Optimum practical and theoretical configurations, circular pins 
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3, Conclusions 

The objectives of this investigation were met by simulating several pin-fin heat 
exchanger configurations with various pin height ratios with an experimentally validated 
three-dimensional finite element model. Specifically: 

1. Comprehensive power function correlations were derived to predict Nusselt 
number and friction factor as a function of Reynolds number, X/D, S/D and H/D for 
3,000<Reo^ <50,000. 

2. Variations in H/D had significant effects on heat exchanger performance. When 
basing performance on wetted area, reducing H/D was found to be most beneficial due to 
the increase in area density. However, when basing performance on heat exchanger face 
area or footprint, increasing H/D provided considerable improvements over the baseline 
configurations especially for streamwise and spanwise pin spacing ratios less than five 
times pin diameter. 

3. By combining all factors found to be beneficial in this study, the optimum 
theoretical configuration was X/D = S/D = 5.0, H/D = 0.25, and X= 0.75 mm. Measuring 
performance on an /? vs. £ chart, this configuration was shown to match the heat transfer 
rate of the baseline case with 300 times less specific energy loss. 

4. When considering practical limitations imposed by particulate contamination 
susceptibility, the best configuration tested was X/D = S/D = 5.0, H/D = 0.25, and X= 6 
mm. This configuration had heat transfer coefficients ranging from 220 to 1050 W/m K 
over a Reynolds number ranging from 50 to 2000 and resulted in over 40 times less 
specific energy loss when compared to the baseline case. 
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D, AIRFOIL-SHAPED PINS 


In this section, the heat transfer and pressure drop characteristics of airfoil-shaped 
pin-fin arrays were investigated using a numerical model. Having thoroughly examined 
the circular pin configurations, changing pin shape was the next logical step in achieving 
performance improvement. 

1, Numerical Simulation 
a. Modeling 

The model simulated a staggered pin array with 10 rows of airfoil-shaped 
pins, but infinitely wide in spanwise extent. As with the circular pin model, this infinite 
span condition was numerically achieved by placing “symmetry”, i.e. zero flux 
boundaries on both the spanwise edges of the array as shown in Figure 44. 

Flow, T = 300K 



boundary, Uy = 0 

Figure 44 Infinite width numerical model, airfoil-shaped pins 

b. Pin Construction Detail 

The pins were geometrically defined by the intersection of two circular 
cylinders as shown in Figure 45. By varying the radius of each circle and controlling 
center-to center distance, several airfoil geometries were obtained. While these two 
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specifications were sufficient to fully define the airfoil shapes, additional parameters 
were used to better characterize the pins as mounted in the test array. These four 
parameters were spanwise ratio {S/D), pin height ratio (H/D), length ratio (LR) and aspect 
ratio (AR). Figure 46 provides detail on how these parameters are defined. For 
simplicity, Dy is replaced by D for spanwise and pin height ratios. 



Figure 45 Airfoil construction 



Aspect Ratio, AR = D^/Dy Length Ratio, LR= D^/X 


Figure 46 Airfoil-shaped pin dimensions 
c. Array Characteristic Length 

The hydraulic diameter definition used for the circular pin phase of study 
was also used for the airfoil-shaped pin-fin array. This was considered to be a useful 
characteristic length as it would account for differences due to variations in H/D, LR and 
AR. 
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d. Grid Independence 

Meshing was carefully specified to ensure maximum coverage in areas of 
interest such as regions of increased velocity and temperature gradient near the walls as 
shown in Figure 47. Additionally, mesh density was meticulously scaled to facilitate grid 
independence studies. The same criteria used to determine grid independence in the 
circular pin study was applied to the airfoil-shaped pin study. Table 12 provides a 
representative example of a grid independence check for a given configuration. 

Table 12 Airfoil array grid independence, S/D = 3, H/D =1, LR =0.7, AR = 2.0, Reoh = 

5,000 


Parameter 

Grid 1 

Grid 2 

Grid 3 

X xy X z 

146/19/11 

158/21/13 

166x22x15 

/pin" 

/13 

/14 

/15 

nodes 

94K 

135K 

163K 

0(W) 

1.51 

1.52 

1.52 

Tou,{^) 

305.37 

305.40 

305.40 

m (kg/s) 

2.81E-04 

2.81E-04 

2.81E-04 













2. Test Approach 

a. Validation 

The numerical model was previously shown to demonstrate satisfactory 
corroboration with historical experiments for circular pins-fins. While the validity of the 
numerical model could be extended to the airfoil-shaped pin-fin array based on this 
corroboration, experiments using actual airfoil-shaped pins were conducted and are 
discussed in a later section. 

b. Test Plan 

A parametric study was devised to investigate the effects of the varying 
S/D, H/D, LR and AR. In this study, each parameter in turn was varied within the ranges 
shown in Table 13. When a single parameter was varied, the other three were maintained 
at nominal values. For each configuration, heat transfer rate and array pressure drop were 
recorded for Reynolds numbers ranging from 3000 to 20000. Based on the trends 
observed, theoretically optimum configurations that combined the best of each parameter 
range were defined and analyzed. 

Table 13 Airfoil numerical test configurations 


Parameter 

Range 

Nominal Value 

S/D 

2-5.0 

3.0 

H/D 

o 

1 

p 

1.0 

LR 

0.5-0.9 

0.7 

AR 

2.4-6.2 

4.2 


3. Results and Discussion 

Upon completion of all the test runs, several key performance indicators were 
examined to understand the heat transfer characteristics and trends for each pin-fin 
configuration and Reynolds number. 

Heat exchanger performance was found to vary significantly as each parameter 
was varied from nominal values. In the sections that follow, Nusselt number and friction 
factor are presented to show the effects of each parameter. Additionally, overall heat 
exchanger performance was evaluated by plotting heat transfer coefficient versus specific 
fluid friction power on a suitable unit surface area basis. 
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a. Effects of Varying Spanwise Spacing 

The effects of varying S/D from 2.0 to 5.0 on heat exchanger performance 
were investigated during the first phase of the analysis. For each configuration, H/D = 
1.0,Ti? = 0.7 and^i? = 4.2. 

Figure 48 provides a plot of Nusselt number as a function of Reynolds 
number for various S/D values. From this chart it is clear that Nusselt number can be 
related to Reynolds number with a power law correlation. Additionally, decreasing S/D 
has a beneficial effect on Nusselt number. However, as shown in Figure 49, friction 
factor also increases as S/D is decreased. Both trends are expected as the decreasing S/D 
provides a more tortuous fluid path and decreased effective flow area. 



Figure 48 Effect of S/D on Nusselt number, airfoil-shaped pins 
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n 



Reynolds number 

Figure 49 Effect of S/D on friction factor, airfoil-shaped pins 


By comparing the configurations on an h vs. E chart as shown in Figure 
50, the overall effects of varying S/D on heat exchanger performance can be determined. 
Here it can be seen that a slight increase in performance can be gained by reducing S/D. 
While the advantage is not overwhelming, it highlights the point that increased heat 
transfer rates associated with smaller S/D outweigh the increase in pressure drop. 

b. Effects of Varying Pin Height Ratio 

Next, the effects of varying H/D from 1.0 to 4.0 on heat exchanger 
performance were studied. For these configurations, S/D = 3.0, LR = 0.7 and AR = 4.2. 

Figure 51 shows a plot of Nusselt number versus Reynolds number for 
various H/D configurations. As in the case of variable S/D, Nusselt number appears to 
follow a power law correlation. However, in this case, the effects of varying H/D are 
more subtle than the effects of variable S/D. Increasing H/D seems to slightly improve 
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Nusselt number. The benefits however would likely fall within the bounds of 
experimental uneertainty and are not significant. 



Friction power/wetted area (W/m^) 

Figure 50 Effects of S/D on h vs. E (wetted area basis), airfoil-shaped pins 

On the other hand, friction factor has a much stronger dependence on pin 
height ratio. As shown in Figure 52, friction factor increases as H/D decreases. This 
trend is mainly due to the geometric differences in the test array. While pressure drop 
decreased as H/D increased, the square of the required average flow velocity to achieve a 
given Reynolds numbers decreased more rapidly resulting in higher friction factors. This 
decrease in velocity was due to the increase in hydraulic diameter as H/D was increased. 
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Nusselt number 



Reynolds number 


Figure 51 Effects of H/D on Nusselt number, airfoil-shaped pins 


n 



3 4 

10 10 

Reynolds number 

Figure 52 Effects of H/D on friction factor, airfoil-shaped pins 
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In order to evaluate the impaet of these trends on overall heat exehanger 
performanee, an h vs. E ehart is shown in Figure 53. When based on wetted area, 
performanee is seen to improve as H/D is deereased. However, when based on faee area 
as shown in Figure 54, an opposite trend is observed. In this ease, performanee is better 
for higher pin ratios. This enhaneement is mainly due to the inereased surfaee area 
provided by the longer pins and the redueed impaet on pin heat transfer coeffieient due to 
endwall effeets. The greater heat rate must be transferred through the heat exehanger 
endwall faee. In theory, longer pins would provide even more performanee. However, 
due to manufaeturing and operating eonstraints, pin height ratios greater than 4.0 are not 
eonsidered representative of eooling sehemes employed in turbine blades. 



10 " 10 ° 10 ^ 10 ^ 


Friction power/wetted area (W/m^) 

Figure 53 Effeet of H/D on h vs. E (wetted area basis), airfoil-shaped pins 
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Figure 54 Effect of H/D on h vs. E (face area basis), airfoil-shaped pins 

c. Effects of Varying Pin-fin Length Ratio 

The effects of varying LR from 0.5 to 0.9 on heat exchanger performance 
were examined. For each configuration, S/D = 3.0, H/D =1.0 and AR = 4.2. 

Figure 55 and Figure 56 show plots of Nusselt number and friction factor 
versus Reynolds number for various pin length ratios. Showing consistent and orderly 
trends, both Nusselt number and friction factor increase as LR is decreased. As before, 
these effects are mainly due to changes in average flow velocity required to achieve 
desired Reynolds number. As length ratio decreases, hydraulic diameter also decreases 
resulting in higher average velocity. Also, pressure drop increases more rapidly than the 
square of the velocity. Thus, the friction factor tends to increase. 
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Figure 55 Effects of length ratio on Nusselt number, airfoil-shaped pins 



Figure 56 Effect of length ratio on friction factor, airfoil-shaped pins 
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The overall effect of LR on heat exchanger performance is realized by 
examining the h vs. E plot shown in Figure 57. From this chart it is evident that the 
increase in heat transfer rate outweighs the increases in pressure drop as LR is reduced. 
Thus reducing pin length ratio improves heat exchanger performance. 



Figure 57 Effect of length ratio on h vs. E (wetted area basis), airfoil-shaped pins 

d. Effects of Varying Pin-fin Aspect Ratio 

Finally, the effects of varying aspect ratio from 2.4 to 6.2 on the 
performance of the heat exchanger were investigated. For these configurations, S/D = 
3.0,/7/D = 1.0 and T7? = 0.7. 

Figure 58 shows a plot of Nusselt number versus Reynolds number for 
each AR variation tested. Interestingly, there is no significant effect on Nusselt number 
as AR is changed. In these cases, the array averaged heat transfer coefficient increases as 
AR increases but is offset by hydraulic diameter which decreases at the same rate. 

Friction factor, on the other hand, is heavily dependent on AR. From Figure 59, it is clear 
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Figure 58 Effects of aspect ratio on Nusselt number, airfoil-shaped pins 



Figure 59 Effects of aspect ratio on friction factor, airfoil-shaped pins 















that friction factor decreases as aspect ratio increases. As increased aspect ratio gives the 
pins a sleeker appearance relative to the flow direetion, this reduetion in friction factor is 
expected. As aspeet ratio is inereased, flow separation is reduced and overall pressure 
drop is dramatically decreased. 

Figure 60 provides a plot of h vs. E for the variable AR ratio runs. It 
eonfirms the expectation that the sleeker or higher aspeet ratio pins provide improved 
performance when compared to low aspeet ratio pins. 



Figure 60 Effect of aspect ratio on h vs. E (wetted area basis), airfoil-shaped pins 

In summary, the four basie defining parameters had a signifieant effect on 

heat exehanger performance. Decreasing S/D, H/D, LR and inereasing AR were all 

beneficial with respect to performance based on wetted area. For performance based on 

face area, the only differenee was that inereasing H/D was most benefieial. Interestingly, 

these trends can be linked to changes in area density, that is, surface area per unit volume. 

Figure 61 shows the effect on area density as eaeh of the four parameters was 
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independently varied from the baseline eonfiguration. It ean be observed that trends in 
heat exehanger performance directly follow area density. Note that area density, like 
performance, increases as S/D, H/D and LR are reduced. Conversely, area density and 
heat exchanger performance increase as AR increases. This linkage can provide helpful 
clues to the heat exchanger designer for preliminary development phases. 



Figure 61 Area density trends for airfoil-shaped pin-fin arrays 

Taking advantage of observed performance trends for each of the four 
defining parameters, it is possible to predict the configuration of an optimum heat 
exchanger within the limits of the test matrix. Such a heat exchanger would have S/D = 
2.0, LR = 0.5 and AR = 6.2. Depending on whether performance was based on wetted 
area or face area would dictate the selection of H/D. For wetted area basis, minimizing 
H/D would be most beneficial whereas for face area basis, maximum H/D would be 
selected. 
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e. Size Effects 

Taking the development one step further, the optimum eonfiguration ean 
be further improved by taking advantage of seale. Figure 62 and Figure 63 show that 
Nusselt number and frietion faetor are essentially independent of array size. This is 
important for two reasons. First, as hydraulie diameter is redueed, heat transfer 
eoeffieient must inerease to maintain Nusselt number. Thus heat transfer rates ean be 
enhaneed by deereasing the axial piteh of an array. Seeondly, the results obtained when 
testing maeroseale heat exehanger arrays ean be direetly applied to mieroseale designs. 
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Figure 62 Effeet of axial piteh on Nusselt number, airfoil-shaped pins 
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Figure 63 Effect of axial pitch on friction factor, airfoil-shaped pins 

f. Optimization 

While size reductions improve heat transfer rate, heat exchanger 

dimensions must include practical limits. Typically, turbine blade designers use a 

minimum passage size of 8 - 12 mils to prevent clogging from airborne particles. Using 

a minimum clearance of 0.3 mm one feasible configuration was S/D = 2, H/D = \, LR = 

0.5 and AR = 6.2 with an axial pitch of 3.75 mm. Area density for this configuration was 

71.5 cm /cm . For this configuration, both passage height and interpin clearance was 0.3 

mm. Another possible configuration would be S/D = 2, H/D = 0.5, LR = 0.5 and AR = 

6.2 with an axial pitch of 7.5 mm yielding an area density of 63 cm /cm . To improve 

performance on a face area basis, the H/D could be raised to 4.0 for either case. 

Numerical simulations indicated that both configurations had nearly identical 

performance. Figure 64 shows an h vs. E plot for an optimum configuration based on 

wetted area. For the sake of comparison, a theoretical worst case, S/D = 5, H/D = 4,LR = 

0.9 and AR = 2.4 is plotted to demonstrate the impact of a poor design. Finally, the 
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optimum and worst-case circular pin configurations from the previous section are 
included to illustrate the advantage gained by switching to the airfoil-shaped pins. 



Friction power/wetted area (W/m^) 


Figure 64 Optimum configuration comparison, airfoil-shaped pins vs. circular pins 

Note that in both the best and worst cases, the airfoil-shaped pins provide 
greater performance advantage as Reynolds number is increased. This is likely due to the 
increased importance of the airfoil shape at higher Reynolds number. Figure 65 shows 
numerically obtained streamline plots for both the airfoil and circularly shaped pins. Both 
pins were in the second row of an array for operating at a Reynolds number of 20,000. 
The airfoil shape tends to minimize boundary layer separation thereby reducing pressure 
drag when compared to circular pins. Since the pressure drag on the pins increases 
proportionally to the square of the flow velocity, the difference in friction power 
associated with the circular pins grows exponentially. 
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Figure 65 Numerical pin-fin streamline comparison, Reoh = 20,000 

In addition to reducing pressure drag, the delayed boundary layer 
separation on the airfoil-shaped pins also provides a heat transfer advantage. Figure 66 
shows a numerically generated comparison of normalized heat transfer coefficient on the 
surface of an airfoil-shaped pin and a circular pin. Note that the heat transfer coefficient 
drops off rapidly on the circular pin as the flow approaches 90 degrees on either side of 
pin. In contrast, the heat transfer coefficient on the airfoil-shaped pin maintains a higher 
value for an additional 30 - 40 degrees. Note also that the normalized heat transfer 
coefficient on the leading edge of the airfoil is significantly greater than it is for the 
circular pin. The actual value in this case was nearly 22 times the average heat transfer 
coefficient on the pin. This was likely due to the sharp leading edge. In this region, the 
boundary layer is extremely thin while the flow velocity is very high. These factors 
combine to yield extremely high local heat transfer coefficients. The graph was not 
scaled to include this portion in order to ensure adequate resolution for characterization 
of the separation point. 
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Figure 66 Local heat transfer coefficient comparison 

The advantage gained by switching to the airfoil shape can be realized by 
comparing the frictional losses between the two optimum configurations of each pin 
shape for a desired heat transfer coefficient. Table 14 displays the comparison between 
the optimum airfoil-shaped pin array and the optimal circular pin-fin array that was 
defined as S/D = 5, X/D = 5, H/D = 0.25 and X = 6 mm. With an area density of 65.6 
cm /cm , the circular pin-fin array closely matched the scale of the airfoil-shaped pin 
array. The second column in the table shows the heat transfer coefficient produced by 
both configurations. Finally the specific fluid friction power required by each 
configuration is compared. In all cases, the airfoil-shaped pin array requires less power 
to produce the desired heat transfer coefficient. Additionally, the relative gains afforded 
by the airfoil design increase at higher Reynolds number. Over the range tested, the 
advantage of switching to the airfoil shapes varies from 15 to 71%. At the highest 
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Reynolds number used in this eomparison, the eireular pin-fin array requires more than 
three times the fiuid power to produee the same heat transfer eoeffieient as the airfoil¬ 
shaped pin array. 


Table 14 Performanee eomparison between eireular and airfoil-shaped pin-fin arrays 


RcDh (airfoil) 

h (W/m^K) 

Ecircular (W/m ) 

£'a/r/oi7(WW) 

% decrease in E 

330 

220 

0.082 

0.070 

15 

900 

450 

1.45 

0.085 

41 

1500 

700 

7.79 

3.97 

49 

2700 

1050 

42.8 

16.9 

60 

4650 

1600 

241 

70.7 

71 


4. Conclusions 

The objeetives of this investigation were met by simulating various airfoil-shaped 
pin-fin heat exehanger eonfigurations with an experimentally validated three-dimensional 
finite element model. The performanee of the airfoil-shaped pin-fin array was found to 
be superior when eompared to similarly sealed eireular pin-fin arrays. This performanee 
inerease was mainly due to the streamlined nature of the airfoil shape whieh was shown 
to delay flow separation from the pins in the array. This not only reduced pressure drag, 
but also increased effective heat transfer surface area on the airfoil-shaped pins. 

Four defining parameters of the airfoil-shaped pins, span wise width ratio, pin 
height ratio, pin length ratio and aspect ratio were varied with a numerical model and the 
effects on heat transfer and pressure drop were explored. The effects of these variations 
tended to follow basic geometric principles. That is, parameter changes that increased 
the area density of an array enhanced overall performance based on heat transfer 
coefficient and fluid friction power when compared to those changes that reduced area 
density. Specifically: 

1. Decreasing S/D was found to increase both Nusselt number and friction 
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factor over the entire Reynolds number range tested. Although frietion faetor was 
inereased, overall performanee was enhanced by redueing span wise pin ratio. 

2. Increasing H/D was found to have little effeet on Nusselt number but 
signifieantly inereased frietion faetor. The inerease in frietion faetor was due to the 
geometry change as pin height ratio inereased. That is, increasing H/D caused hydraulic 
diameter to increase significantly thereby reducing average velocity requirement for a 
given Reynolds number. Thus friction factor increased. The overall impact of changes in 
H/D depended on whether wetted area or face area was used as a performance basis. 
Reducing H/D improved performanee on a wetted area basis. Conversely, increasing 
H/D enhanced performanee based on faee area whieh may be important for eooling 
eleetronie eomponents. 

3. As in the ease of pin height ratio, changes in pin length ratio signifieantly 
affeeted hydraulic diameter and thus affeeted Nusselt number and friction factor. 
Reducing length ratio increased both Nusselt number and frietion factor. Overall 
performanee was also enhaneed when length ratio was redueed. 

4. Pin aspect ratio was the only parameter that had an opposite effeet on area 
density and overall performance. When aspeet ratio was inereased, frietion faetor was 
redueed as the pins beeame more streamlined. No signifieant ehanges oeeurred in 
Nusselt number so overall performanee was improved as aspeet ratio was inereased. 

5. By eombining the faetors shown to improve performance into one design, 
an optimum eonfiguration was developed. Sized to resist partieulate eontamination, the 
optimum eonfiguration was S/D =2, H/D = 0.5, LR = 0.5, AR = 6.3 and X = 6 mm. 
Compared to a circular pin configuration with approximately the same area density, the 
airfoil-shaped pin-fin array mueh lower frictional losses to maintain a given heat transfer 
coefficient. The advantage of the airfoil-shaped pin array grew from 15 to 71% as 
Reynolds number was inereased from 330 to 4650. 
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III. EXPERIMENTAL ANALYSIS 


A, OVERVIEW 

This section covers the experimental testing completed to support the numerical 
analysis presented in the previous sections. Up to this point, the numerical model had 
only been validated against historical data from Metzger et al. (Ref 5) and Arora (Ref. 
11). Based on these comparisons, the model was found to be in excellent agreement with 
experimental heat transfer data. Friction factor predictions were not as accurate but were 
considered useful for determining trends between various pin-fin array configurations. 
No experimental data were available to validate the airfoil-shaped pins. In this section, 
experimental testing of three circular and one airfoil-shaped pin-fin array will be 
discussed. 

B, EXPERIMENTAL SETUP 

1. Experimental Apparatus 

The experimental rig consisted of a rectangular entrance duct, pin-fin test array, 
mixing chamber and exit duct. Airflow was driven by a two-stage regenerative blower 
mounted at the end of the exit duct. A schematic of the test apparatus layout is shown in 
Figure 67. 

The entrance duct was constructed of plexiglass and measured three meters in 
length to ensure a fully developed velocity profile at the test section. The inner 
dimensions of the duct were 33 mm x 250 mm to match the test section. The duct was 
open to ambient air and was fitted with a screen to prevent debris from entering the test 
section and fouling the downstream blower. All duct joints were carefully constructed to 
maintain airtight integrity and present a smooth surface to the incoming air thereby 
preventing flow disturbances. 

The test section was comprised of 10 separate sections of pins. Figure 68 shows 
test section construction detail. Each section or row consisted of aluminum pin-fins 
mounted with screws between two aluminum plates material measuring 12 mm x 50 mm 
X 250 mm.. Each section was thermally insulated from the adjacent section by 1 mm 
plexiglass strips to facilitate independent row-by-row heat transfer analysis. The 
adiabatic sidewalls were made of 12 mm thick plexiglass and permitted excellent 
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Figure 67 Experimental apparatus sehematie 

visibility into the test array. As with the entranee duet, eare was taken to minimize flow 
disturbanees eaused by uneven endwall alignment. Additionally, rubber gasket material 
was used in all joints to prevent air leaks in the test array. Using a smoke generator, all 
leaks were revealed and subsequently eorreeted before testing was eondueted. 
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Figure 68 Test seetion detailed diagram 

For some of the testing, two different meehanieal turbulenee generators were 
installed immediately prior to the test seetion. One consisted of wide metallic strips and 
the other was a standard screen mesh. Figure 69 shows photographs of these generators. 

Flexible electric heating elements were bonded to the endwalls using a thermal 
adhesive. Rated at 50W each, these heaters were designed to meet the heating 
requirements for the entire Reynolds number range to be tested. A simple program was 
written to provide automatic on/off control of these heaters to maintain desired endwall 
temperature and to record total heat rate based on electrical energy requirements. As a 
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Figure 69 Turbulence generators, coarse (left) and fine (right) 

backup means of measuring total heat flux, an in-line wattmeter was used to provide an 
additional source for measuring electrical power supplied to the strip heaters. Wall 
temperature of each section was indicated by thermocouples mounted in grooves 
machined into the endwalls. For this application, the thick aluminum plate served to 
provide adequate thermal inertia to promote an even temperature distribution spatially 
and temporally throughout the endwall. 

Other instrumentation measured air temperature, pressure and velocity. Air 
temperature was recorded using a single thermocouple at the entry duct inlet and four 
thermocouples at the end of the insulated mixing chamber. Air pressure was recorded at 
the test section inlet and exit. Airflow was measured using a turbine flowmeter. Flow 
velocity was controlled using relief valves mounted between the blower and the airflow 
meter and by placing filters on the inlet. Further details of the experimental configuration 
can be found in Ref. 20, and photographs of some of the equipment are shown in 
Appendix D. 
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2, Experimental Method 

Test procedures were consistently and methodically executed to ensure reliable 
data quality. At the start of each test period, ambient atmospheric pressure was recorded 
using NOAA data from a nearby airport located at the same altitude as the laboratory. 
Next, the strip heaters were energized and set to maintain desired wall temperature. 
Before turning on the blower, heater energy was measured to calculate heat losses to the 
environment. Finally the blower was activated and various bleed valve and entry filter 
combinations were used to achieve desired flow rates. Steady state conditions were 
achieved when heat exchanger outlet temperatures remained constant. Streaming data 
were recorded in a spreadsheet format and included individual wall temperatures, cooling 
air inlet temp, outlet air temperature and pressure, turbine flowmeter output and pressure 
drop across the array. Data were recorded for 20 minute periods to ensure sufficient data 
sampling. 

a. Array Characteristic Length 

For consistency, volume-based hydraulic diameter was used for the 
experimental testing as it had been for the numerical analysis. 

b. Test Matrix 

Several configurations of pin-fin arrays were chosen for testing as shown 
in Table 15. Both experimental and numerical runs were completed with each 
configuration. These configurations were chosen to provide a sample of data points 
applicable to each of the major numerical analysis categories previously completed. 
Additionally, pin dimensions were selected to ensure compatibility with the common test 
section and ducting to reduce setup time and cost. Each configuration was tested with a 6 
and 12 degree C temperature difference between wall temperature and cooling air 
temperature and had an axial pitch of 50 mm. Finally, each configuration was also tested 
with and without turbulence generators installed immediately upstream of the test section 
to qualitatively characterize the effects of turbulence intensity on heat exchanger 
performance. 
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Table 15 Experimental and numerieal test eonfigurations, Reoh = 5,000 - 45,000 


Circular Pin Configurations 

Diameter 

S/D 

X/D 

H/D 

33.0 mm 

1.5 

1.5 

1.0 

16.5 mm 

3.0 

3.0 

2.0 

10.0 mm 

5.0 

5.0 

3.33 

Airfoil Pin Configuration 


S/D 

H/D 

LR 

AR 

3.67 

2.44 

0.9 

3.30 


3, Data Reduction 

Once the raw data were acquired, a MATLAB program was used to calculate 
desired output quantities. The following steps were followed to calculate the test array 
Nusselt number and friction factor. 

a. Nusselt Number 

Total heat transfer rate and mass flow rate were recorded directly from the 
instrumentation and data acquisition system. 

Outlet bulk temperature was calculated from a simple energy balance as: 


T = r. + - 

out in 


Q 


m C„ 


( 26 ) 


As a check, this outlet temperature was compared with the outlet bulk 
temperature obtained by averaging the results of the four thermocouples located in the 
mixing chamber. In most cases, the temperatures were identical. If the temperatures 
differed by more than 0.1 K, then the run was repeated. With output temperature 
resolved, the log mean temperature difference was calculated over the test array as 
follows: 


_ {T„a„ -Tbulk,„) - - Tbulk ^^,) 

KaU-Tbulk, '^ 

The heat transfer coefficient was calculated using the following 


h 


Finally, the Nusselt number for the array was determined using: 
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Nu^ 

k 

Additionally, due to the independent nature of the modular heated 
seetions, row-by-row temperature profde was obtained through a simple energy balanee. 
From this, row-by-row Nusselt number development was derived. 
b. Friction Factor 

The frietion faetor in the array was ealeulated using the following 
AP D, 

definition; f= ^ 

\ pU^L 

where L was the overall streamwise length of the array. 

4. Uncertainty Analysis 

Experimental uneertainty was estimated using the method of Kline and 
MeClintoek (Ref 21) for Nusselt number, frietion faetor and Reynolds number. Due to 
the wide range of Reynolds numbers tested, the relative uneertainty values also eover a 
wide range. Table 16 shows these estimated values for eaeh parameter. For Reynolds 
number and frietion faetor, the greatest uneertainty oeeurred at the lowest Reynolds 
number and the smallest uneertainty oeeurred at the high Reynolds numbers. Of 
partieular note is the high uneertainty values for the airfoil-shaped pins. This was 
primarily due to the diffieulty assoeiated with measuring the extremely low pressure drop 
aeross the array at the lowest Reynolds number tested even with an inelined manometer. 
Details on the uneertainty analysis are provided in Appendix E. 

Table 16 Experimental Uneertainty 


Parameter 

Relative Uneertainty (%) 

(5,000-45,000) 

4.9- 1.9 


8.3 

NuDh {AT= 12K) 

4.2 

Frietion faetor (33 mm pins) 

3.3-0.62 

Frietion faotor(16.5 mm pins) 

3.4-0.64 

Frietion faotor(10 mm pins) 

20.0-0.68 

Frietion faetor (airfoil pins) 

25 - 0.78 
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C. EXPERIMENTAL RESULTS 

Upon completion of the test runs, the experimental results were eompared with 
data obtained from the numerieal model. As diseussed previously, the numerieal models 
were used to develop Nusselt number and frietion faetor eorrelations as a funetion of 
Reynolds number. Additionally, heat transfer eoeffieient was plotted against fluid 
frietion power as a key performanee metrie. In this seetion, results from eaeh 
experimental eonfiguration are matehed with the eorresponding numerieal results to 
provide a means for evaluating the model validity. 

1, Nusselt Number 

Figure 70 shows Nusselt number results for eaeh of the pin eonfigurations. Note 
that in most oases there is exoellent agreement between the numerieal and experimental 
results for the mid to low range Reynolds numbers. However, for Reynolds numbers 
above 20,000, the numerieal model tends to result in higher Nusselt number values when 
eompared to the experimental data. The largest differenoe noted was 18%. There are at 
least two possible reasons for this differenoe. 

The first oause may be due to oontaot resistanoe between the endwall and pin-fins 
in the experimental rig. In the numerieal model, both the endwall and pin-fin surfaoes 
were maintained at the same temperature for all runs. However, in the experimental rig, 
the pin-fins are heated indireotly by the strip heaters via oonduotion through the endwalls. 
At the lower Reynolds numbers, the pin-fins are able to maintain desired temperature as 
the heat flux requirement is relatively small. However, at the higher Reynolds numbers, 
the small yet finite oontaot resistanoe results in a slightly lower pin temperature when 
eompared to the endwall. This differenoe in temperature oan result in lower Nusselt 
numbers than predioted by the numerieal model. 

The seoond reason may be due to inaoouraoies in the numerieal solution at higher 
Reynolds numbers. Although turbulenoe models are generally more aoourate at high 
Reynolds number, they often have diffioulties oharaoterizing regions of flow reoiroulation 
suoh as those that would form behind the pin-fins in the test array. Additionally, there 
may be problems in aoourately depioting points of flow separation from the pins. It is 
reasonable to eonsider that inaoouraoies arising from poor oharaoterization of these 

regions would be magnified as flow velooity inoreased. 
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However, despite the differenees noted, the agreement between the numerieal 
model and the experimental results were considered satisfactory and within the bounds of 
experimental uncertainty. 



10 ^ 10 '* 10 ^* 10 ^ 10 '* 10 ^ 
ReDh ReDh 
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ReDh 


Figure 70 Model validation, Nusselt number 

Due to the ability to heat each row independently, it is possible to resolve row-by- 
row local Nusselt number. Figure 71 and Figure 72 show local Nusselt number plotted 
for each row of pins. Both figures show the 33 mm pin configuration. In both cases, the 
Nusselt number quickly climbs to maximum value by the third row. However, Nusselt 
number decreases rapidly in Figure 71 similar to the heat transfer coefficient trends seen 
in the numerical case represented by Figure 23. However, in Figure 72, the Nusselt 
number remains high after the third row. This may be due to the increased turbulent 
energy present in the higher speed flow. The numerical model also showed this type of 
level Nusselt development curve at high flow velocities. These similarities increased 
confidence in the numerical model’s ability to characterize the heat transfer 
characteristics of the pin-fin heat exchanger. 
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2, Friction Factor 

Figure 73 shows the eomparison between frietion faetor for the experimental and 
numerieal runs for eaeh eonfiguration. For the 33 mm diameter eireular pin-fin array, the 
numerieal model overestimates the frietion faetor by up to 50%. In the ease of the 
smaller eireular pins and airfoil-shaped pins, the agreement is mueh better. However, in 
all eases exeept for the 10 mm pins, the numerieal solution and experimental results 
appear to diverge as Reynolds number inereases. This differenee at high Reynolds 
number is eonsistent with the trends observed above with Nusselt number. Unlike the 
heat transfer solution however, only the aforementioned eoneem with the numerieal 
modeling of the reeireulation zone is relevant. 
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Figure 73 Model validation, frietion faetor 

In faet, possible diffieulties in numerieally eharaeterizing the reeireulation zones 

or flow separation points are reasonable eonsidering the trends observed with frietion 

faetor. The size of the reeireulation zones relative to the open volume would be greater 
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for the larger pins than for the small eireular or airfoil-shaped pins. Additionally, these 
zones would have a greater impaet at higher Reynolds number for all eonfigurations. 
These same trends were observed in the earlier comparison with the experimental results 
of Arora (Ref. 11). 

Further evidence of premature flow separation can be seen in the numerical 
solution. In the previous section, a numerically derived plot of normalized heat transfer 
coefficient development on both a circular and airfoil-shaped pin-fin were shown as 
Figure 66. Here, the heat transfer coefficient appeared to drop off at 80 - 85 degrees from 
the leading edge of the circular pin which is indicative of flow separation. This 
separation point is significantly earlier than the typically accepted value for turbulent 
flow of approximately 110 - 120 degrees as found in elementary fluid mechanics 
textbooks such as White (Ref. 22). As a result of the earlier separation, the low pressure 
recirculation zone or wake is likely larger for the numerical solution than it is for the 
experimental case. This would cause higher pressure drop and subsequently higher 
friction factors for the circular pins. 

For the airfoil-shaped pins, there is very little separation and the recirculation 
zone is not as significant as it was for the circular pin case. As this zone is relatively 
small, its contribution to the total pressure drop in the array is also small. As a result, 
errors associated with the numerical solution are also not very significant. As a result, 
there is significantly better corroboration with the experimental solution. 

Despite the problems with accurately predicting the pressure drop for the large 
diameter circular pin array, the numerical model was still considered useful in showing 
trends between different configurations. In Figure 74, two plots of friction factors are 
shown for all configurations. One chart includes only numerical results and the other 
only shows experimental results. Although the numerical model overestimates the large 
pin friction factor, it can be seen that the experimental trends are properly predicted. 
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Figure 74 Friction factor trends between configurations 

















3, Performance Comparisons 

Plots of heat transfer eoeffieient versus speeifie fluid frietion power provide a 
very useful metrie in determining the performanee of heat exehanger. These eharts 
illustrate the “eost” of a given eonfiguration to transfer heat. Figure 75 shows these 
performanee eharts for eaeh eonfiguration using data from the numerieal and 
experimental runs. Note that despite the previously identified differenees in frietion 
faetor seen with the eireular pins, the h vs. E eurves show good agreement between 
numerieal and experimental results. This is likely due to the faet that the model tends to 
overestimate both Nusselt number and pressure drop at the higher Reynolds numbers. 
These errors are effeetively eaneeled sinee the higher heat transfer eoeffieient appears to 
“eost” more in terms of fluid frietion power. This eauses the numerieal performanee 
eurve to shift slightly up and to the right eompared to the experimental data for the same 
Reynolds number range. These slight differenees did not impair the model’s utility in 
eomparing the relative performanee between eonfigurations. 
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Figure 75 Model validation, overall performanee 
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4, Effect of Turbulence Generators 

A qualitative evaluation of the effect of turbulence generators was conducted and 
the results are shown in this section. Figure 76 shows the effects of the turbulence 
generators on each pin-fin type. For the circular pins, there was no noticeable difference 
in Nusselt number when the turbulence generators were installed. However, the 
turbulence generators increased Nusselt number in the airfoil-shaped pin-fin array by 35 
- 50%. This difference is likely due to the difference in flow characteristics between 
these configurations. With the circular pins, there is already a high level of turbulence 
within the array once the flow gets past the first or second row. Therefore introducing 
additional turbulence upstream of the array has little effect on the already highly 
turbulent flow situation. However, in the case of the airfoil-shaped pins, the flow 
naturally has a lower turbulence level as the path is not as tortuous as it is in the circular 
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Figure 76 Effect of turbulence generators on Nusselt number 
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pin case. As a result, introdueing upstream disturbanees signifieantly inereases 
turbulenee intensity within the array yielding an enhaneed heat transfer rate. It is 
interesting to note that there was little differenee in Nusselt number between the eoarse 
and fine turbulenee generators for the airfoil-shaped pins. This suggests a turbulenee 
threshold beyond whieh there is little benefit in ereating larger flow disturbanees. 

While the introduetion of the turbulenee generators affeeted the experimental 
results, no similarities were observed in the numerieal model. Inereasing turbulenee 
intensity in the numerieal model had no signifieant effeet on the heat transfer rate or 
pressure drop. This suggests one of two possibilities. First, like the eireular pin 
experimental ease, the baekground turbulenee in the numerieal airfoil-shaped pin model 
may have been suffieiently high sueh that inereased turbulenee intensity had little or no 
effeet. Second, the level of turbulent intensity achieved numerieally may have been 
mueh lower and not representative of the level attained by the turbulenee generators in 
the experiment. In either ease, the previous seetion detailing the numerieal analysis of the 
airfoil-shaped pins eovers a nominal situation in whieh the flow, although turbulent, is 
relatively undisturbed as it enters the array. In a representative applieation like turbine 
blade eooling, the flow would likely be much more disturbed. Based on the experimental 
impaet of the turbulenee generators, aetual heat transfer rates for installed airfoil-shaped 
pin-fin arrays may be even greater than predieted by the results of this numerieal work. 

Unfortunately, due to the eonfiguration of this experiment, no meaningful data 
eould be obtained to evaluate the effeets of turbulence generators on frietion faetor. This 
was beeause the turbulenee generators were loeated between the entry pressure tap and 
the test array. Thus reeorded pressure drop would have ineluded the turbulenee 
generators rather than just the test array. 

D, CONCLUSIONS 

The objeetives of this phase of testing were met by eomparing heat transfer and 
pressure drop results obtained by the numerieal model with data reeorded from various 
pin configurations mounted in an experimental rig. Several findings were: 
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1. The numerical simulation modeled the heat transfer characteristics of the 
circular and airfoil pin-fin arrays over a Reynolds number range of 5,000 to 45,000 to 
within 18%. 

2. The numerical model overestimated friction factor by up to 50% for the 33 
mm circular pin array but was much closer for the 16.5 mm and 10 mm circular pin and 
airfoil-shaped pin arrays. These differences may have been due to shortcomings in the 
simulation’s ability to properly model flow separation and areas of recirculation. In most 
cases, the numerically obtained friction factor tended to diverge from the experimental 
results at higher Reynolds number. 

3. Despite shortcomings in characterizing absolute friction factor, the 
numerical model was useful for discriminating between configurations since relative 
friction factor values were consistent. 

4. Overall heat exchange performance based on heat transfer coefficient and 
fluid friction power was characterized closely by the numerical model for all 
configurations. 

5. The introduction of turbulence generators in the experimental rig 
increased Nusselt number by 35-50% for the airfoil-shaped pin-fin array. No significant 
changes were noted for the circular pin arrays. This was likely due to the fact that the 
circular pin arrays were already operating at high turbulence intensity due to the tortuous 
nature of the flow around the pin-fins. 

6. The numerical model was not affected by the introduction of additional 
turbulence. This suggests that the airfoil-shaped pin array may perform even better than 
predicted in a real-world environment. 
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IV. CONCLUSIONS 


Due to the modular nature of this study, specific conclusions have already been 
discussed at the end of each major section/chapter. However, this section shall be used to 
review the significant findings of this work as an entity and capture the main concepts in 
summary form. 

The major contributions of this study were to advance knowledge in the growing 
field of microscale heat exchanger technology and to develop innovative ways of 
improving heat transfer performance in the staggered cross pin-fin heat exchanger. The 
major objective of this work was met by successfully simulating the pin-fin heat 
exchanger with a 3-D numerical model and using it to investigate a wide range of 
configurations, Reynolds numbers, array sizes, pin height ratios and pin shapes. The 
model was validated with historical and current experimental data. It consistently 
predicted Nusselt number for all configurations tested within 18% but was found to have 
some shortcomings in predicting friction factor. Despite these shortcomings, it provided 
a reliable measure of overall heat exchanger performance. 

The model demonstrated the advantages of reducing axial pitch to achieve 
superior heat exchanger performance within the limits of particulate contamination 
resistance. Using a porous medium type open flow volume based hydraulic diameter as 
the characteristic length, the model demonstrated that Nusselt number and friction factor 
were independent of actual heat exchanger size for all configurations, pin height ratios 
and pin shapes. This important finding permits microscale heat exchanger optimization 
to be performed on a more manageable macroscale basis since the dimensionless 
correlations can be directly transferred to the small scale. 

Additionally, the model demonstrated that airfoil-shaped pin-fin arrays 
outperform similar sized circular pin-fin arrays. In the optimum airfoil-shaped pin array 
case modeled, the circular pin-fin array required three times the specific energy loss to 
produce the same heat transfer rate. In addition, experimental data have indicated that 
increased turbulence levels can improve heat transfer rates in airfoil pin arrays by 35- 
50% while having no benefit for the circular pin array. Thus the performance of the 
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airfoil-shaped pin array, already superior to the cireular pin array, would be further 
enhaneed in an installed eonfiguration where turbulenee was likely to be greater than in 
the laboratory or numerieal environment. 

Further, optimum pin configuration depends greatly on the desired application. In 
situations requiring large heat flux through a small face area, such as in electronic 
component cooling applications, increases in pin height ratio have a significant effect on 
performance. This is due to the magnification of wetted by pin-fin arrays thereby 
increasing effective heat transfer coefficient. Details on this area magnification and a 
sample application are presented in Appendix F. The trends observed in this study would 
be useful to provide the designer with initial guidance in producing the optimum heat 
exchanger design for a given application. 

Finally, it is important to understand some limits to the applicability of the 
correlations presented in this work. First, since the effects of buoyancy and body forces 
were not investigated, results in a rotating reference frame may vary. As such, direct use 
of these correlations would be more suited to stator blades or fixed electronic component 
cooling applications. However, even for rotor blades, these correlations would provide 
the design engineer with a means of comparing the relative benefits of various pin-fin 
configurations. Secondly, even if concerns over particulate contamination could be 
eliminated, there is a limit in size reduction based on the mean free path of the gas. 
Beyond this limit, the gas would no longer obey the continuum hypothesis on which the 
Navier-Stokes equations are based and could violate the no-slip condition along rigid 
boundaries. Appendix G provides further details on this limitation. 
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APPENDIX A. MODEL DETAILS 

GENERAL METHOD 

All models were eonstrueted by first ereating a three-dimensional “bloek” to 
represent the eooling air in the test array. The pins were ereated by removing 
appropriately shaped volumes using Boolean subtraetion from the initial bloek. This 
resulted in holes in the model that represented the pin-fins in the test array. Finally two 
additional 3-D bloeks representing the adiabatie entry and exit duet were ereated and 
“glued” to the test array. 

ROUND PIN MODELS 

The round pins were eonstrueted in ANSYS by simply ereating eireular volumes 
with the desired radius. 

AIRFOIL-SHAPED PINS 

The airfoil-shaped pins were ereated from the interseetion of two eireular volumes 
of radius, R, as shown in Figure 77. The horizontal eomponent of radius, Rx, was 
defined as a fraetion of axial piteh, X, by speeifying the desired length ratio (2*Rx/X) 
from the test matrix. Then, slenderness ratio (Ry/Rx) was used to define Ry and achieve 
the desired aspect ratio (A/D). Finally, R was calculated and the circle centers were 
spaced 2Ry apart. The volume defined by the intersection of these circles was subtracted 
from the cooling air block and represented the airfoil-shaped pins. 



Figure 77 Airfoil construction details 
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ISOTHERMAL PIN BOUNDARY CONDITION JUSTIFICATION 

As stated earlier, the pins were treated as isothermal entities through boundary 
eondition seleetion. Caleulations to justify this deeision are present in this seetion. The 
eonfiguration used for these ealeulations is: S/D = X/D = 1.5, H/D = 1, X= 50 mm. This 
was one of the experimental eonfigurations. 


The governing equations for this oaleulation are taken from Ineropera and DeWitt 
(Ref 18). Using an adiabatie tip assumption, they are: 


0 T-T^ eoshm(T-x) 
T^-T^ cosh (ml) 


[hP 

where m = I - , is the cooling air temperature, T is the pin base 

temperature, and T is the temperature along the pin at a distance x from the base. Also, L 
is half the pin height (0.033 m) in order to use the adiabatic assumption. 

Perimeter, P = 0.1046 m, kaiuminum =177 (W/mK), 

and cross sectional area, Ac = 8.7092e-4 

For high Reynolds number (where this analysis is particularly important) h can be 
taken as approximately 150 W/m K, based on experimental results. 

Substituting, m = 8.32 


Thus, r = 7;+(7;-7;) 


cosh 8.32(0.0165-x) 
cosh(8.32* 0.0165) 


Assuming worst case conditions, = 300A and 7), = 312A 

The midpin temperature is 311.89K. This is considered close enough to 312K to 
assume that the pin is isothermal. 
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APPENDIX B. SAMPLE ANSYS MACRO 


In order to efficiently and systematically manage the large number of numerical 
models evaluated in this study macros were written for each heat exchanger type. Each 
macro performed the following basic tasks: 


1. Provided easy input entry such as Reynolds number, pin spacing and 
height ratio, axial pitch, inlet air and wall temperature, reference pressure and mesh 
density. 

2. Constructed and meshed the model. 

3. Set ANSYS solution parameters to include solver type and iteration 
number. 


4. Ran the solution. 


In reviewing the following sample macro, the reader will notice that the solution 
is started three times. This is to accommodate the sensitive nature of the Collocated 
Galerkin (COLG) solver. For the initial solution, density is held constant and the MSU 
approach is used to discretize the momentum equations. No thermal analysis is 
conducted. This allows a rough flow solution to develop. For the next series of iterations, 
density is allowed to vary according to the equation of state and the energy equation is 
solved. For the third phase, COFG is used to discretize only the momentum equation. 
By this point, the solution is fairly robust and the final phase can be initiated. In the final 
phase, the COFG approach is used for all parameters. In some cases, if the preliminary 
steps are not completed, the solution will diverge. The sample macro is presented next: 
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/TITLE,Re=5000,X/D=1.5,S/D=1.5Jw =306,h/v/p =7.8/15/2,SKE,COLG 
adv,10atm,1.5inf 

! creates an endless 1.5 pin fin array with 10 rows 
! unheated entry and exit regions 
! L.J. Hamilton 7 NOV 02 

I****************************************** 

lEnter Reference pressure in Pa 
Pref= 101350 

lEnter your Reynold's Number (VanFossen definition) 

Re=5000 

lEnter your X/D (cylinder spacing in the flow direction) 
xd=1.5 

lEnter your S/D (spanwise (perp to flow) 
sd=1.5 

lEnter your H/D (Pin height/diameter ratio) 
hd=1.0 

lEnter xy (h)/and z (v) grid scaling 
h = 7.8 
V = 15 

lEnter pin factor 

p = 2 

lEnter axial pitch (Absolute distance between pin centers) 
x=.0127 

lEnter your Twall 
Twall=306 

Tfilm = (Twall+300)/2 


I calculates pin diameter 
d=x/xd 

ICalculate kinematic viscosity using Sutherland Law 
rho = Pref/(287*Tfilm) 
rho300 = Pref/(287*300) 

Imu = 1.86e-5 

mu = 383.4/(Tfilm+110.4)*1.71e-005*(Tfilm/273)**(3/2) 
nu = mu/rho 

Icalculated for a unit cell (s*2x) 
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Varray = 2*sd*xd*hd*d**3 
Vpin = 3.14159*hd*d**3/4 
Vopen = Varray - 2*Vpin 
Awpin = 3.14159*hd*d**2 
Afeet = 3.14159*d**2/4 
Awall = 2*sd*xd*d**2 - 2*Afeet 
Aw = 2*Awpin + 2*Awall 
Abar = Vopen/(2*x) 

Aduct = sd*hd*d**2 
Dprime = 4*Vopen/Aw 


i*******************|^l0^ Velocity Csiculstion********* 

Vin = Re*mu*Abar/(Dprime*rho300*Aduct) 

,*******************Pg^gi.^.^g length of pin section***** 

!Pin center to pin center streamwise distance 
ddx=d*xd 

!Pin center to pin center spanwise distance 
ssy=d*sd 
!Pin Height 
Ht = hd*d 

ICaicuiate entry iength for FD turb fiow 

EntryDh = 2*Ht 

EntryRe = Vin*EntryDh/nu 

entry = EntryDh*4.4*EntryRe**(1/6) 

entryinit = entry 

lEnter exit iength (after the pins. Metz:7.62cm, to 
! prevent feedback in the ANSYS 

! soln: 12.5cm)[m] 

exit = 0.125 

ITotai iength of pin section with (xd/2)*d before the 
lieading edge and (xd/2)*d foiiowing the traiiing edge 
!of the iast row of pins. 
xiength= (ddx*10) 

ITotai width of ANSYS model, insulated wall to sym plane. 
ylength= (ssy*1.5) 

|*********gg^y ******************************** 

, ****************jl^lg i^oDULES SETS FLOTRAN rarmS********* 
!***********selects operating preference*************** 

/NOPR 
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/PMETH,OFF,0 

KEYW,PR_SET,1 

KEYW,PR_STRUC,0 

KEYW,PR_THERM,1 

KEYW,PR_FLUID,0 

KEYW,PR_ELMAG,0 

KEYW,MAGNOD,0 

KEYW,MAGEDG,0 

KEYW,MAGHFE,0 

KEYW,MAGELC,0 

KEYW,PR_MULTI,1 

KEYW,PR_CFD,1 

/GO 

!* 

/COM, 

/COM,Preferences for GUI filtering have been set to display: 
/COM, Thermal 
/COM, FLOTRANCFD 

!* 

/UNITS,SI 

/PREP7 

!* 

Iselects element type 
ET,1,FLUID142 

!* 

i*****************l^ljiy^l_ 50LN OPTIONS ************* 

FLDATA1,SOLU,TRAN,0 

FLDATA1,SOLU,FLOW,1 

FLDATA1,SOLU,TEMP,0 

FLDATA1,SOLU,TURB,1 

FLDATA1,SOLU,COMP,0 

FLDATA1,SOLU,VOF,0 

FLDATA1 ,SOLU,SFTS,0 

FLDATA1,SOLU,IVSH,0 

FLDATA1,SOLU,SWRL,0 

FLDATA1 ,SOLU,SPEC,0 

FLDATA1,SOLU,ALE,0 

,******************!EXECUTION CONTROL ******** 
/COM„Steady State Analysis,0 
FLDATA2,ITER,EXEC,30, 

FLDATA2,ITER,OVER,0, 

FLDATA2,ITER,APPE,0, 

FLDATA3,TERM,VX,0.01, 

FLDATA3,TERM,VY,0.01, 

FLDATA3,TERM,VZ,0.01, 
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FLDATA3,TERM,PRES,1e-008, 

FLDATA3,TERM,TEMP,1e-008, 

FLDATA3,TERM,ENKE,0.01, 

FLDATA3,TERM,ENDS,0.01, 

FLDATA5,OUTP,SUMF,10, 

FLDATA5,0UTP,YPLU,T 

,******************^PP^P MAKE DENSITY AN OUTPUT 
FLDATA5,OUTP,DENS,T, 

!* 

/PREP7 

FINISH 

/PREP7 

!************!NITIAL FLUID PROPERTIES (CONSTANT DENSITY)************* 

FLDATA12,PROP,DENS,4 

FLDATA13,VARY,DENS,0 

FLDATA12,PROP,VISC,4 

FLDATA13,VARY,VISC,1 

FLDATA12,PROP,COND,4 

FLDATA13,VARY,COND,1 

FLDATA12,PROP,SPHT,4 

FLDATA13,VARY,SPHT,1 

!* 

FLDATA7,PROT,DENS,AIR-SI 

FLDATA8,NOMI,DENS,-1 

FLDATA9,COF1,DENS,0 

FLDATA10,COF2,DENS,0 

FLDATA11 ,COF3,DENS,0 

FLDATA7,PROT,VISC, AIR-SI 

FLDATA8,NOMI,VISC,-1 

FLDATA9,COF1,VISC,0 

FLDATA10,COF2,VISC,0 

FLDATA11 ,COF3,VISC,0 

FLDATA12,PROP,IVIS 

FLDATA7,PROT,COND,AIR-SI 

FLDATA8,NOMI,COND,-1 

FLDATA9,COF1,COND,0 

FLDATA10,COF2,COND,0 

FLDATA11 ,COF3,COND,0 

FLDATA7,PROT,SPHT, AIR-SI 

FLDATA8,NOMI,SPHT,-1 

FLDATA9,COF1,SPHT,0 

FLDATA10,COF2,SPHT,0 

FLDATA11 ,COF3,SPHT,0 

!***Select ref temp = 300 K************* 

FLDATA14,TEMP,NOMI,300 


105 



FLDATA14JEMPJTOT,300 
FLDATA14,TEMP,BULK,300 
FLDATA15,PRES,REFE,Pref 

!* 

!***PLACEHOLDER FOR TURBULENCE MODEL INPUTS************** 


!SKE=1(Default),ZETM=2,RNG=3,NKE=4,GIR=5,SZL=6 

!FLDATA24,TURB,MODL,1 

!*************TEMP CFD SOLVER IS PGMR ******************** 

FLDATA18,METH,TEMP,4 

FLDATA22,MAXI,TEMP,1000, 

FLDATA20,SRCH,TEMP,12, 

FLDATA20A,PGMR,FILL,6, 

FLDATA20A,PGMR,MODP,0, 

FLDATA21,CONV,TEMP,1e-12, 

FLDATA23,DELT,TEMP,1e-010, 

!* 


!***********ADVECTION PARMS 


*************** 


FLDATA,ADVM,MOME,MSU 

FLDATA,ADVM,TURB,MSU 

FLDATA,ADVM,PRES,MSU 

FLDATA,ADVM,TEMP,MSU 

!* 


I************************ BLOCK & PINS ********************* 


/prep7 

ICreates initial block with extra 0.01 length 


iStart in middle of first row 
xloc=entry + ddx/2 
start=entry 
*do,i,1,5 

Block,entry,entry+ddx,0,ylength,0,Ht/2 

cyl4,xloc,ssy/2,d/2„„Ht/2 

cyl4,xloc,(ssy*1.5),d/2„„Ht/2 

Block,entry+ddx,entry+2*ddx,0,ylength,0,Ht/2 

cyl4,(xloc+ddx),0,d/2„„Ht/2 

cyl4,(xloc+ddx),ssy,d/2„„Ht/2 

xloc=xloc + (ddx*2) 

entry=entry + (ddx*2) 


*enddo 


I**************** SUBTRACT PINS ******************** 

FLST,2,10,6,ORDE,10 

FITEM,2,1 

FITEM,2,4 
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FITEM,2,7 

FITEM,2,10 

FITEM,2,13 

FITEM,2,16 

FITEM,2,19 

FITEM,2,22 

FITEM,2,25 

FITEM,2,28 

FLST,3,20,6,ORDE,20 

FITEM,3,2 

FITEM,3,-3 

FITEM,3,5 

FITEM,3,-6 

FITEM,3,8 

FITEM,3,-9 

FITEM,3,11 

FITEM,3,-12 

FITEM,3,14 

FITEM,3,-15 

FITEM,3,17 

FITEM,3,-18 

FITEM,3,20 

FITEM,3,-21 

FITEM,3,23 

FITEM,3,-24 

FITEM,3,26 

FITEM,3,-27 

FITEM,3,29 

FITEM,3,-30 

VSBV,P51X,P51X 

!************ADD ENTRY BLOCK ******************* 
Block,0,start,0,ylength,0,Ht/2 

!************ADD EXIT BLOCK ******************** 


Block,start+xlength,start +xlength+exit,0,ylength,0,Ht/2 
Ineed to glue entry and exit! 

FLST,2,12,6,ORDE,4 

FITEM,2,1 

FITEM,2,-2 

FITEM,2,31 

FITEM,2,-40 

VGLUE,P51X 

!***Model is now built 


*******************jjp,gj^gg jggg^g^^jg******** 
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N\E\N, 1 ,1,1,1 
/ANG, 1 
/REP, FAST 
/AUTO, 1 
/REP 
VPLOT 

/color,pbak,off 

/REPLOT 


I************ PINS ****************************** 

!1/n 1/n 1/n 1/n 1/n 1/n 1/n 1/n 1/n 1/n 1/n 1/n 

!***********Cylinders (1/n) 
lsel,s,line„16 

*do,i,0,1 

*do,k,0,4 

iCylinders exit/upper side 
lsel,a,line,,16 + 1*i + 64*k 
iCylinders inlet/lower side 
lsel,a,line,,56 + 5*i + 64*k 
*enddo 
*enddo 


*do,i,0,1 

*do,j,0,1 

*do,k,0,4 

iCylinders exit/upper side 

lsel,a,line„48 + 1*i + 10*j + 64*k 

ICylinders inlet/lower wall side 

lsel,a,line„14 + 5*i + 10*j + 64*k 

*enddo 

*enddo 

*enddo 


LESIZE,all,„3*h/xd,1/p„„1 


Innnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnn 


*********************** 


Cylinders (n) 


lsel,s,line„15 

*do,i,0,1 

*do,k,0,4 

ICylinders inlet/upper side 
lsel,a,line,,15 + 3*i + 64*k 
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ICylinders exit/lower wall side 
lsel,a,llne,,55 + 7*1 + 64*k 
*encldo 
*enddo 

*do,1,0,1 
*do,j,0,1 
*do,k,0,4 

ICylinders Inlet/upper side 

lsel,a,llne„47 + 3*1 + 10*j + 64*k 

ICylinders exit/lower wall side 

lsel,a,llne„13 + 7*1 + 10*j + 64*k 

*enddo 

*enddo 

*enddo 

LESIZE,all,„3*h/xd,1/p„„1 

l************Test Section Symmetry lines, short**************** 

lsel,s,llne„321 

lsel,a,llne„322 

lsel,a,llne„219 

lsel,a,llne„220 

lsel,a,llne„228 

lsel,a,llne„229 

lsel,a,llne„281 

lsel,a,llne„282 

lsel,a,llne„301 

lsel,a,llne„302 

lsel,a,llne„307 

lsel,a,llne„308 

lsel,a,llne„365 

lsel,a,llne„366 

lsel,a,llne„194 

lsel,a,llne„199 

lsel,a,llne„173 

lsel,a,llne„174 

*do,l,0,3 

lsel,a,llne„153 + I 
lsel,a,llne„327 + I 
lsel,a,llne„335 + I 
lsel,a,llne„343 + I 
lsel,a,llne„351 + I 
*enddo 
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LESIZE,all„,h*xd/3.1 


,******************jgg^ Section Symmetry Lines, Long*************** 

isei,s,line„130 

isei,a,line„135 

isei,a,line„164 

isei,a,line„165 

isei,a,iine„217 

lsei,a,line„218 

isei,a,line„237 

isei,a,iine„238 

isei,a,line„243 

isei,a,line„244 

lsei,a,line„283 

isei,a,line„284 

isei,a,iine„292 

isei,a,iine„293 

isei,a,line„361 

isei,a,line„362 

isei,a,line„363 

isei,a,line„364 

isei,a,iine„179 

isei,a,line„180 

LESIZE,ali,„3*h.1 


I**************** 


Entry/Exit Symmetry Lines******************** 


isei,s,line„7 

isei,a,line„36 

isei,a,iine„90 

isei,a,line„91 


LESIZE,aii,„15*h,2 


9 9 9 9 


1 


!********Entry/Exit Symmetry Lines for proper spacing 

isei,s,line„71 

isei,a,line „100 

isei,a,iine„26 

isei,a,line„27 


********** 


LESIZE,ali,„15*h,1/2„„ 


1 


1 **************!^l 0 ^ snd Exit Linos******************* 


no 






^ | 0 ^ SHCl ExIt LlflGS******************* 

lsel,s,line„2 

lsel,a,line„37 

lsel,a,line„89 

lsel,a,line„92 
LESIZE,all„,3*h*sd/xd.1 

l***************Q I U© ^ j|-|0g************************* 

lsel,s,line„25 

lsel,a,line„28 

lsel,a,line,,3 

lsel,a,line,,6 

lsel,a,line„35 

lsel,a,line„38 

lsel,a,line„67 

lsel,a,line„70 

lsel,a,line„99 

lsel,a,line„102 

lsel,a,line„131 

lsel,a,line„134 

lsel,a,line„163 

lsel,a,line„166 

lsel,a,line„195 

lsel,a,line„198 

lsel,a,line„227 

lsel,a,line„230 

lsel,a,line„259 

lsel,a,line„262 

lsel,a,line„66 

lsel,a,line„101 

LESIZE,all„,3*h*sd/xd.1 

!*****Vertical Lines Entry + those needed to ensure proper mesh*********** 

lsel,s,line„52 

lsel,a,line„45 

lsel,a,line„109 

lsel,a,line„116 

LESIZE,all„,v,1/5„„1 

i***********Yg|.^jjjgl Lines, Test Section + Exit*** 

lsei,s,line„110 

isei,a,line„115 

isei,a,iine„51 


111 






Isel,a,line„46 


*do,i,0,1 

*do,k,0,8 

lsel,a,line,,10 + 1*i + 32*k 

*enddo 

*enddo 


,*****************^gi^.^gl (pins)* 

*do,i,0,1 

*do,j,0,1 

*do,k,0,4 

lsel,a,line„21 + 1*i + 10*j + 64*k 

lsel,a,line„53 + 1*i + 10*j + 64*k 

*encldo 

*enddo 

*enddo 


LESIZE,all„,v,5„„1 


lsel,all 

asel,all 


**************** 


***************** 

asel,s,area„9 


Boundary Conditions 
Pins vx=vy=vz=0, T=Twaii 


************* 


******* 


*do,k,0,4 

*do,i,0,1 

asel,a,area„9 + i + 28*k 
asei,a,area„27 + i + 28*k 
*enddo 

asei,a,area„14 + 28*k 
asei,a,area„23 + 28*k 
*enddo 


DA,aii,vx,0,1 
DA,ail,vy,0,1 
DA,all,vz,0,1 
DA,all,TEMP,Twall,1 
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I********************* 


Entry/Exit endwalls vx=vy=vz=HFLU=0********** 


asel,s,area,,1 

asel,a,area„12 


DA,all,vx,0,1 
DA,all,vy,0,1 
DA,all,vz,0,1 
SFA,ALL„HFLUX,0 


I ********************** 11 


, vx=vy=vz=0, T=Twall***** 


asel,s,area,,25 

asel,a,area„32 

asel,a,area„49 

asel,a,area,,54 

asel,a,area„69 

asel,a,area,,77 

asel,a,area„63 

asel,a,area„82 

asel,a,area,,91 

asel,a,area„41 

DA,all,vx,0,1 

DA,all,vy,0,1 

DA,all,vz,0,1 

DA,all,TEMP,Twall,1 


!*******symmetry (side) vy=0, HFLU=0 

asel,s,area,,7 

asel,a,area„24 

asel,a,area„141 

asel,a,area,,30 

asel,a,area,,45 

asel,a,area„149 

asel,a,area,,53 

asel,a,area,,60 

asel,a,area,,157 

asel,a,area„68 

asel,a,area„73 

asel,a,area,,165 

asel,a,area„81 

asel,a,area„88 

asel,a,area„173 

asel,a,area„39 

asel,a,area,,16 

asel,a,area,,15 

asel,a,area„36 
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asel,a,area,,40 

asel,a,area„86 

asel,a,area„170 

asel,a,area„80 

asel,a,area„72 

asel,a,area„162 

asel,a,area„67 

asel,a,area,,58 

asel,a,area,,154 

asel,a,area„52 

asel,a,area,,44 

asel,a,area„146 

asel,a,area„29 

asel,a,area„22 

asel,a,area,,4 

DA,all,vy,0,1 
SFA,ALL„HFLUX,0 

!**********symmetry (middle) vz=0, hfLU=0****** 

asel,s,area„2 

asel,a,area„26 

asel,a,area,,35 

asel,a,area,,50 

asel,a,area,,57 

asel,a,area,,64 

asel,a,area,,71 

asel,a,area,,78 

asel,a,area„85 

asel,a,area,,92 

asel,a,area„43 

asel,a,area,,13 

DA,all,vz,0,1 

SFA,ALL„HFLUX,0 

!************!nlet vx=Vin vy=vz=0, t=300K********* 

asel,s,area,,8 

DA,all,vx,Vin,1 

DA,all,vy,0,1 

DA,all,vz,0,1 

DA,all,TEMP,300,0 


PClOf~0 

asel,s,area„21 
DA,all,pres,0,0 


*********** 
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asel,all 

I **************** jY|0g j-|***************************** 

FLST,5,12,6,ORDE,2 

FITEM,5,1 

FITEM,5,-12 

CM,_Y,VOLU 

VSEL, , , ,P51X 

CM,_Y1,V0LU 

CHKMSH,'VOLU' 

CMSEL,S,_Y 

!* 

VSWEEP,_Y1 

!* 

CMDELE,_Y 

CMDELE,_Y1 

CMDELE,_Y2 

!* 

/Ul,MESH,OFF 

!***********run FIRST (0-30 iter) solution*************** 

/SOLU 

FINISH 

/SOLU 

SOLVE 

!*********SECOND SOLN (30-60 iter) OPTIONS************* 

FLDATA1,SOLU,TRAN,0 

FLDATA1,SOLU,FLOW,1 

FLDATA1,SOLU,TEMP,1 

FLDATA1,SOLU,TURB,1 

FLDATA1,SOLU,COMP,0 

FLDATA1,SOLU,VOF,0 

FLDATA1 ,SOLU,SFTS,0 

FLDATA1,SOLU,IVSH,0 

FLDATA1,SOLU,SWRL,0 

FLDATA1 ,SOLU,SPEC,0 

FLDATA1,SOLU,ALE,0 

!* 

!*********second execution controls************* 

/COM„Steady State Analysis,0 
FLDATA2,ITER,EXEC,30, 

FLDATA2,ITER,OVER,0, 

FLDATA2,ITER,APPE,0, 

FLDATA3,TERM,VX,0.01, 

FLDATA3,TERM,VY,0.01, 
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FLDATA3,TERM,VZ,0.01, 

FLDATA3,TERM,PRES,1e-08, 
FLDATA3,TERM,TEMP,1e-08, 
FLDATA3,TERM,ENKE,0.01, 

FLDATA3,TERM,ENDS,0.01, 

FLDATA5,OUTP,SUMF,10, 

FLDATA5,OUTP,YPLU,T 

,******************^PP^P jQ MAKE DENSITY AN OUTPUT 
FLDATA5,OUTP,DENS,T, 

!* 

I************ SECOND FLUID PROPERTIES *************** 

FLDATA12,PROP,DENS,4 

FLDATA13,VARY,DENS,1 

FLDATA12,PROP,VISC,4 

FLDATA13,VARY,VISC,1 

FLDATA12,PROP,COND,4 

FLDATA13,VARY,COND,1 

FLDATA12,PROP,SPHT,4 

FLDATA13,VARY,SPHT,1 

!* 

FLDATA7,PROT,DENS,AIR-SI 

FLDATA8,NOMI,DENS,-1 

FLDATA9,COF1,DENS,0 

FLDATA10,COF2,DENS,0 

FLDATA11 ,COF3,DENS,0 

FLDATA7,PROT,VISC, AIR-SI 

FLDATA8,NOMI,VISC,-1 

FLDATA9,COF1,VISC,0 

FLDATA10,COF2,VISC,0 

FLDATA11 ,COF3,VISC,0 

FLDATA12,PROP,IVIS 

FLDATA7,PROT,COND,AIR-SI 

FLDATA8,NOMI,COND,-1 

FLDATA9,COF1,COND,0 

FLDATA10,COF2,COND,0 

FLDATA11 ,COF3,COND,0 

FLDATA7,PROT,SPHT, AIR-SI 

FLDATA8,NOMI,SPHT,-1 

FLDATA9,COF1,SPHT,0 

FLDATA10,COF2,SPHT,0 

FLDATA11 ,COF3,SPHT,0 

!* 

I*********** SECOND ADVECTION PARMS *************** 

FLDATA,ADVM,MOME,MSU 

FLDATA,ADVM,TURB,MSU 
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FLDATA,ADVM,PRES,MSU 

FLDATA,ADVM,TEMP,MSU 

!* 

!********run SECOND (30-60 iter) solution************* 

/SOLU 

FINISH 

/SOLU 

SOLVE 

!*******BEGIN THIRD solution (60-90 iter) options****************** 
!*********THIRD SOLN OPTIONS ************* 

FLDATA1,SOLU,TRAN,0 

FLDATA1,SOLU,FLOW,1 

FLDATA1,SOLU,TEMP,1 

FLDATA1,SOLU,TURB,1 

FLDATA1,SOLU,COMP,0 

FLDATA1,SOLU,VOF,0 

FLDATA1 ,SOLU,SFTS,0 

FLDATA1,SOLU,IVSH,0 

FLDATA1,SOLU,SWRL,0 

FLDATA1 ,SOLU,SPEC,0 

FLDATA1,SOLU,ALE,0 

!* 

i*********jl^ll^Q EXECUTION CONTROLS************* 
/COM„Steady State Analysis,0 
FLDATA2,ITER,EXEC,30, 

FLDATA2,ITER,OVER,0, 

FLDATA2,ITER,APPE,0, 

FLDATA3,TERM,VX,0.01, 

FLDATA3,TERM,VY,0.01, 

FLDATA3,TERM,VZ,0.01, 

FLDATA3,TERM,PRES,1e-08, 

FLDATA3,TERM,TEMP,1e-08, 

FLDATA3,TERM,ENKE,0.01, 

FLDATA3,TERM,ENDS,0.01, 

FLDATA5,OUTP,SUMF,10, 

FLDATA5,OUTP,YPLU,T 

,******************^PP^P MAKE DENSITY AN OUTPUT 
FLDATA5,OUTP,DENS,T, 

!* 

!************THIRD FLUID PROPERTIES ************* 

FLDATA12,PROP,DENS,4 

FLDATA13,VARY,DENS,1 

FLDATA12,PROP,VISC,4 

FLDATA13,VARY,VISC,1 

FLDATA12,PROP,COND,4 
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FLDATA13,VARY,C0ND,1 

FLDATA12,PROP,SPHT,4 

FLDATA13,VARY,SPHT,1 

!* 

FLDATA7,PR0T,DENS,AIR-SI 

FLDATA8,N0MI,DENS,-1 

FLDATA9,COF1,DENS,0 

FLDATA10,COF2,DENS,0 

FLDATA11 ,COF3,DENS,0 

FLDATA7,PR0T, Vise, AIR-SI 

FLDATA8,N0MI,VISC,-1 

FLDATA9,COF1,VISC,0 

FLDATA10,COF2,VISC,0 

FLDATA11 ,COF3,VISC,0 

FLDATA12,PR0P,IVIS 

FLDATA7,PR0T,C0ND,AIR-SI 

FLDATA8,N0MI,C0ND,-1 

FLDATA9,COF1,COND,0 

FLDATA10,COF2,COND,0 

FLDATA11 ,COF3,COND,0 

FLDATA7,PR0T,SPHT, AIR-SI 

FLDATA8,N0MI,SPHT,-1 

FLDATA9,COF1,SPHT,0 

FLDATA10,COF2,SPHT,0 

FLDATA11 ,COF3,SPHT,0 

!* 

!***********THIRD ADVECTION PARMS *************** 

FLDATA,ADVM,MOME,COLG 

FLDATA,ADVM,TURB,MSU 

FLDATA,ADVM,PRES,MSU 

FLDATA,ADVM,TEMP,MSU 

!* 

!*******run THIRD (60-90 iter) solution*********** 

/SOLD 

FINISH 

/SOLD 

SOLVE 

!************FINAL (90 - max iter) solution****** 
!*********FINAL SOLN OPTIONS ************* 

FLDATA1,SOLU,TRAN,0 

FLDATA1,SOLU,FLOW,1 

FLDATA1,SOLU,TEMP,1 

FLDATA1,SOLU,TURB,1 

FLDATA1,SOLU,COMP,0 

FLDATA1,SOLU,VOF,0 
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FLDATA1 ,SOLU,SFTS,0 

FLDATA1,SOLU,IVSH,0 

FLDATA1,SOLU,SWRL,0 

FLDATA1 ,SOLU,SPEC,0 

FLDATA1,SOLU,ALE,0 

!* 

!*********FINAL EXECUTION CONTROLS ************* 

/COM„Steady State Analysis,0 
FLDATA2,ITER,EXEC,120, 

FLDATA2,ITER,OVER,0, 

FLDATA2,ITER,APPE,0, 

FLDATA3,TERM,VX,0.01, 

FLDATA3,TERM,VY,0.01, 

FLDATA3,TERM,VZ,0.01, 

FLDATA3,TERM,PRES,1e-08, 

FLDATA3,TERM,TEMP,1e-08, 

FLDATA3,TERM,ENKE,0.01, 

FLDATA3,TERM,ENDS,0.01, 

FLDATA5,OUTP,SUMF,10, 

FLDATA5,OUTP,YPLU,T 

,******************^PP^P jQ MAKE DENSITY AN OUTPUT 
FLDATA5,OUTP,DENS,T, 

!* 

!************FINAL FLUID PROPERTIES ********************************* 

FLDATA12,PROP,DENS,4 

FLDATA13,VARY,DENS,1 

FLDATA12,PROP,VISC,4 

FLDATA13,VARY,VISC,1 

FLDATA12,PROP,COND,4 

FLDATA13,VARY,COND,1 

FLDATA12,PROP,SPHT,4 

FLDATA13,VARY,SPHT,1 

!* 


FLDATA7,PROT,DENS,AIR-SI 
FLDATA8,NOMI,DENS,-1 
FLDATA9,COF1,DENS,0 
FLDATA10,COF2,DENS,0 
FLDATA11 ,COF3,DENS,0 
FLDATA7,PROT,VISC, AIR-SI 
FLDATA8,NOMI,VISC,-1 
FLDATA9,COF1,VISC,0 
FLDATA10,COF2,VISC,0 
FLDATA11 ,COF3,VISC,0 
FLDATA12,PROP,IVIS 
FLDATA7,PROT,COND,AIR-SI 
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FLDATA8,N0MI,C0ND,-1 
FLDATA9,COF1,COND,0 
FLDATA10,COF2,COND,0 
FLDATA11 ,COF3,COND,0 
FLDATA7,PR0T,SPHT, air-si 
FLDATA8,N0MI,SPHT,-1 
FLDATA9,COF1,SPHT,0 
FLDATA10,COF2,SPHT,0 
FLDATA11 ,COF3,SPHT,0 
!* 

I*********** FINAL ADVECTION FARMS *************** 

FLDATA,ADVM,M0ME,C0LG 

FLDATA,ADVM,TURB,C0LG 

FLDATA,ADVM,PRES,C0LG 

FLDATA,ADVM,TEMP,C0LG 

!* 


i**********************i-u^ FINAL (90“ msx itGr) 


solution 


********************************** 


/SOLU 

FINISH 

/SOLU 

SOLVE 


I ********************* QO to IdSt 


g0^************************ 


/POST1 
FINISH 
/POST1 
SET, LAST 


/post1 
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APPENDIX C. DISCUSSION OF ANSYS SOLVERS 


The text in this appendix was taken from the ANSYS Theory reference manual 
(Ref 17) to provide the reader with deeper insight into the algorithms used by 
FLOTRAN. The equation numbers in this section refer to equations from that reference 
and are not intended to be in numerical sequence with the equations in the body of the 
dissertation. 

DERIVATION OF FLUID FLOW MATRICES 


A segregated, sequential solution algorithm is used. This means that element matrices are 
formed, assembled and the resulting system solved for each degree of freedom separately. 
Development of the matrices proceeds in two parts. In the first, the form of the equations 
is achieved and an approach taken towards evaluating all the terms. Next, the segregated 
solution algorithm is outlined and the element matrices are developed from the equations. 

Discretization of the Equations 

The momentum, energy, species transport, and turbulence equations all have the form of 
a scalar transport equation. There are four types of terms; transient, advection, diffusion, 
and source. For the purposes of describing the discretization methods, let us refer to the 

variable considered as . The form of the scalar transport equation is: 

Equation 7,68, 


_a 

at 

_a 

9k 


^ ^dx) 9yl^ ^9y^ 


9zr ^9z 


+ S. 


(pVjC/) = 


where: 
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c 

r 

s 


'^= transient and advection coefficient 
diffusion coefficient 
'^= source terms 


Transport Equation Representation below shows what the variables, coefficients, and 
source terms are for the transport equations. The pressure equation is derived using the 
continuity equation. Its form will unfold during the discussion of the segregated solver. 
The terms are defined in the previous section. 

Since the approach is the same for each equation, only the generic transport equation 
need be treated. Each of the four types of terms will be outlined in turn. Since the 
complete derivation of the discretization method would require too much space, the 
methods will be outlined and the reader referred to more detailed expositions on the 
subjects. 

Transport Equation Representation 



Meaning 

DOE 

c* 

r* 


Vx 

x-velocity 

VX 

1 


P& - 3 

p/3x + Rx 

Vy 

y-velocity 

VY 

1 

e 

5 

p/^y + Ry 

Vz 

z-velocity 

1 

VZ 

1 


Pgz - 3 

p/3z + Rz 

T 

temperature 

TEMP 

Cp 

K 

Qv + 

+ p. (J) 3 
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p/^t 

k 

kinematic 

energy 

ENKE 

1 



s 

dissipation 

rate 

ENDS 

1 


sis** 

Y. 

species mass 

fraction 

SPOl- 

06 

1 

P 

Dmi 

0 


The discretization process, therefore, consists of deriving the element matrices to put 
together the matrix equation; 

Equation 7,69, 

[[Af+ [aI"''®"''”" 1 + [a"''"'”"]){ 4.e) - {s*} 

Galerkin's method of weighted residuals is used to form the element integrals. Denote by 
W® the weighting function for the element, which is also the shape function. 

Transient Term 

The first of the element matrix contributions is from the transient term. The general form 
is simply: 

Equation 7,70, 
j-^tjansientj ^ 

A lumped mass approximation is used so that 

Equation 7,71, 
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^(^d(vol) = %*lfW= 

at ^ ^ at J 

A backward difference scheme is used to evaluate the transient derivative. On a nodal 
basis, the following implicit formulation is used. The current time step is the nth time 
step and the expression involves the previous two time step results. 

Equation 7.72, 

ajpt) (pkn-? 

at 2At 2At 2At 

For a Volume of Fluid (VOF) analysis, the above equation is modified as only the results 
at one previous time step are needed: 

Equation 7.73. 

^(PP _ (Pt)n (Ptf., 

at At At 

The above first-order time difference scheme is chosen to be consistent with the current 
VOF advection algorithm. 

The nth time step produces a eontribution to the diagonal of the element matrix, while the 
derivatives from the previous time step form contributions to the source term. 

Advection Term 

Currently Flotran has three approaches to discretize the advection term. The monotone 
streamline upwind (MSU) approach is first order accurate and tends to produce smooth 
and monotone solutions. The streamline upwind/Petro-Galerkin (SUPG) approach and 
the collocated Galerkin (COLG) approach are second order accurate and tend to produce 
oscillatory solutions. 
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Monotone Streamline Upwind Approach (MSU) 


The advection term is handled through a monotone streamline approach based on the idea 
that pure advection transport is along characteristic lines. It is useful to think of the 
advection transport formulation in terms of a quantity being transported in a known 
velocity field. See Streamline Upwind Approach. 

Figure 7.1, Streamline Upwind Approach 


Downstream Value 



The velocity field itself can be envisioned as a set of streamlines everywhere tangent to 
the velocity vectors. The advection terms can therefore be expressed in terms of the 
streamline velocities. 

In pure advection transport, one assumes that no transfer occurs across characteristic 
lines, i.e. all transfer occurs along streamlines. Therefore one may assume that the 
advection term. 

Equation 7.74, 

dx dz ds 
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when expressed along a streamline, is eonstant through out an element: 


Equation 7,75. 


r^ad,ecl,onj . jVV®d(vol) 


This formulation is made for every element, each of which will have only one node 
which gets contributions from inside the element. The derivative is calculated using a 
simple difference: 

Equation 7,76, 

d(pC^Vs) _ -(pC/s<l^)p 

ds Ds 


where: 


D = subscript for value at the downstream node 

U = subscript for value taken at the location at which the streamline through the 
downwind node enters the element 

^ s = distance from the upstream point to the downstream node 

The value at the upstream location is unknown but can be expressed in terms of the 
unknown nodal values it is between. See Streamline Upwind Approach again. 

The process consists of cycling through all the elements and identifying the downwind 
nodes. A calculation is made based on the velocities to see where the streamline through 
the downwind node came from. Weighting factors are calculated based on the proximity 
of the upwind location to the neighboring nodes. 

Additional details are provided by Rice and Schnipke (Ref 23). 
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Streamline Upwind/Petro-Galerkin Approach (SUPG) 

The SUPG approach consists of a Galerkin discretization of the advection term and an 
additional diffusion-like perturbation term which acts only in the advection direction. 

Equation 7.77. 

r^advection ] . .yya ^ ^ 

L ^ J J 9k: 9/ 9z 

^ , zh fv.aw® Vy9/V® V,9W® 

C2j( --+ -i-+ ^- 

^^J2U^3g 9k 9/ 9z 

Vx9(pC/) ^ Vy9(pC^(^) ^ V;9(pC^(^) 

9k 9y 9z 

where: 

Co 

global coefficient set to 1.0 
h = element length along advection direction 

Umag = 

1 if 0 < Pe < 3 

Pe/3 ifPe>3 

Pe = ^ = Peclet number 

2r^ 

It is clear from the SUPG approach that as the mesh is refined, the perturbation terms 
goes to zero and the Galerkin formulation approaches second order accuracy. The 
perturbation term provides the necessary stability which is missing in the pure Galerkin 
discretization. Additional details are provided by Brooks and Hughes (Ref 24). 





127 



Collocated Galerkin Approach (COLG) 


The COLG approach uses the same discretization scheme with the SUPG approach with 
a collocated concept. In this scheme, a second set of velocities, namely, the element- 
based nodal velocities are introduced. The element-based nodal velocities are made to 
satisfy the continuity equation, whereas the traditional velocities are made to satisfy the 
momentum equations. 


’^advection' ^ jy^e 


3(pv'C.^<|.) 3(pv'C/) 3(pv|C.j4,) 


dK 




dL 


cl(vol) 




zh 


2U® 




9k 


di 


dz 


j v|a(pc/) v;a(pC4) v|a(pc/) ]^^ 


dK 


dz 


Where all the parameters are defined similar to those in the SUPG approach. 

In this approach, the pressure equation is derived from the element-based nodal 
velocities, and it is generally asymmetric even for incompressible flow problems. The 
collocated Galerkin approach is formulated in such a way that, for steady-state 
incompressible flows, exact conservation is preserved even on coarse meshes upon the 
convergence of the overall system. 

Diffusion Terms 

The expression for the diffusion terms comes from an integration over the problem 
domain after the multiplication by the weighting function. 

Equation 7.79. 


128 



Diffusion contribution = JW 

J 


afp an 




dx 

d' 
dz 


‘^dx 




di 


cl(vol) 


cl(vol) 


\ / 

The X, y and z terms are all treated in similar fashion. Therefore, the illustration is with 
the term in the x direetion. An integration by parts is applied; 

Equation 7.80, 


JW® 


dx 


an 


A ax 




Onee the derivative of ^ is replaced by the nodal values and the derivatives of the 
weighting function, the nodal values will be removed from the integrals 


Equation 7,81. 

f 

ax 


Equation 7,82. 

SW® 


w® = 


dx 


The diffusion matrix may now be expressed as: 

Equation 7,83, 

I'^^diffusionj ^ jw^r^W® +W®r^W® +Wz%W® d(vol) 
Source Terms 


129 



The evaluation of the source terms consists of merely multiplying the source terms as 
depicted in Streamline Upwind Approach by the weighting function and integrating over 
the volume. 

Equation 7,84. 

Sj-JW'S^cl(vol) 

Segregated Solution Algorithm 

Each degree of freedom is solved in sequential fashion. The equations are coupled, so 
that each equation is solved with intermediate values of the other degrees of freedom. 

The process of solving all the equations in turn and then updating the properties is called 
a global iteration. Before showing the entire global iteration structure, it is necessary to 
see how each equation is formed. 

The preceding section outlined the approach for every equation except the pressure 
equation, which comes from the segregated velocity-pressure solution algorithm. In this 
approach, the momentum equation is used to generate an expression for the velocity in 
terms of the pressure gradient. This is used in the continuity equation after it has been 
integrated by parts. This nonlinear solution procedure used in FLOTRAN belongs to a 
general class of Semi-Implicit Method for Pressure Linked Equations (SIMPLE). There 
are currently two segregated solution algorithms available. One is the original SIMPLEF 
algorithm, and the other is the enhanced SIMPLEN algorithm. 

The incompressible algorithm is a special case of the compressible algorithm. The change 
in the product of density and velocity from iteration to the next is approximating by 
considering the changes separately through a linearization process. Denoting by the 
superscript * values from the previous iteration, in the x direction, for example, results: 

Equation 7,85, 
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pVx =pVx + P*Vx -p*Vx 

The continuity equation becomes: 

Equation 7,86, 


dp , . Q(p^«) , . ^(pS) . 

dt dK dx d)/ 

3(p*vz]^3(pvi) 9 (pX) a(p\;) a(p\*)_ 
dz dz dx dz 


0 


The transient term in the continuity equation can be expressed in terms of pressure 
immediately by employing the ideal gas relationship: 

Equation 7,87, 

fW® ■^cl(vol) = — fW® -^cl(vol) 

J 91 91J RT 


The backward differencing process is then applied directly to this term. 
Application of Galerkin's method to the remaining terms yields: 

Equation 7,88, 
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JW 


a(p\,) ^ a(p\) a(p\ 

dK dz 


■JW 


-JW 


a 

(pv*x) ^ 


ax 

; 


a 



ax 

L 


9/ 9z 


cl(vol) 


cl(vol) 




di 


d(vol) 


There are thus three groups of terms. In the first group, terms with the derivatives of the 
unknown new velocities must be integrated by parts to remove the derivative. The 
integration by parts of just these terms becomes: 


Equation 7,89, 


JW 


0(p'Vx) 3(p*Vy) 3(p*Vz) 


Ox 


0/ 


dz 


d(vol) 


= JW ^p*Vx + P*Vy +p*V2j d(area) 


-J 



,aA/ / * 

\ aA/ / * 1 

.dH' 

(P Vx] 

'a< fiP ^ 


' dz 


d(vol) 


Illustrating with the x direction, the unknown densities in the second group expressed in 
terms of the pressures are: 


Equation 7,90, 


Jw|(pv*x)d(vol).J 


W_^ 
R ax 



In the third group, the values from the previous iteration are used to evaluate the 
integrals. 
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The next step is the derivation of an expression for the veloeities in terms of the pressure 
gradient. When the momentum equations are solved, it is with a previous value of 
pressure. Write the algebraie expressions of the momentum equations assuming that the 
eoeffieient matriees eonsist of the transient, adveetion and diffusion eontributions as 
before, and all the souree terms are evaluated exeept the pressure gradient term. 

Equation 7.91. 

E 

AVx=S^-XW ^ Cl(vol) 

e=1 

Equation 7.92. 

E 


Equation 7.93. 



Eaeh of these sets represents a system of N algebraie equations for N unknown veloeities. 
It is possible, after the summation of all the element quantities, to show an expression for 
eaeh veloeity eomponent at eaeh node in terms of the veloeities of its neighbors, the 
souree terms whieh have been evaluated, and the pressure drop. Using the subseript "i" to 
denote the nodal equation, for i = 1 to N, where N is the number of fluid nodes and 
subseript "j" to denote its neighboring node: 

For SIMPLEF algorithm: 

Equation 7.94. 
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Equation 7.95, 



_ 1 _ 

II 



Equation 7.96, 



For SIMPLEN algorithm; 







d(vol) 


cl(vol) 


cl(vol) 


where for SIMPLEF algorithm: 

Equation 7.100, 
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"Z^ijVxj+Sx 


J 



Equation 7.101. 



Equation 7.102. 


j^i 

j_ 


and or SIMPLEN algorithm: 


Vxi - 


_J_ 

a- j’‘i 
-i^ + Va" 

fX ^ u 


^Yi = 


-Zaij(Vyj-Vyi)+biy 

_ 

a-y J*i 
^+va.y 

V ■ U 
J ^ 
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j 


3,1 j’^i 
_!L+y 

j 


dtl 


Here the ay represent the values in the x,y and z coefficient matrices for the three 
momentum equations, r is the relaxation factor, and hi is the modified source term taking 
into effect the relaxation factors. 


For the purposes of this expression, the neighboring velocities for each node are 
considered as being known from the momentum equation solution. At this point, the 
assumption is made that the pressure gradient is constant over the element, allowing it to 
be removed from the integral. This means that only the weighting function is left in the 
integral, allowing a pressure coefficient to be defined in terms of the main diagonal of the 
momentum equations and the integral of the weighting function: 

For SIMPLEF algorithm: 


Equation 7.106. 


^ii e=1 


Equation 7.107. 


My 



N 

^Wcl(vol) 


Equation 7.108. 


. N 

Mj = — cl(vol) 


e=1 
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For SIMPLEN algorithm; 


Mx 



r 


X 


1 


j^i 

E 

j 


a-f 


N 

2]Wd(vol) 

e=1 


My 


1 



j^i 

y a.y 
j 


N 

^Wd(vol) 

e=1 


Mj 


_!L + y a /i- 
j 'J 


j^i 


N 


^Wd(vol) 

e = 1 


Therefore, expressions for unknown nodal veloeities have been obtained in terms of the 
pressure drop and a pressure coefficient. 


Equation 7.112, 


Vx 


= Vx-Mx 




Equation 7.113. 


V 


y 



di 


Equation 7.114, 



These expressions are used to replace the unknown velocities in the continuity equation 

to convert it into a pressure equation. The terms coming from the unknown velocities 

137 



(replaced with the pressure gradient term) and with the unknown density (expressed in 
terms of the pressure) contribute to the coefficient matrix of the pressure equation while 
all the remaining terms will contribute to the forcing function. 


The entire pressure equation can be written on an element basis, replacing the pressure 
gradient by the nodal pressures and the derivatives of the weighting function, putting all 
the pressure terms on the left hand side and the remaining terms on the right hand side 
(equation 7.115). 

Equation 7.115, 


[pf j 


oW 






m sw 


p*M 


m dAI 


dK d)/ ^ di di 




m 


dz 


d(vol)‘ 


i 

T 


W 


A 

’'T 




* Pi a f * Pi 

v„- +— V,- 


T di\ 


T 


d(vol)' 


sw ... aw aw 

P P Vy +—p 


ax 




az 


d(vol)' 


qw 




a/ 


Vz 


d(vol)' 


-Jw[p*Vx] d(area)® - JW ^p*Vy j d(area)® - JW ^p^v^j d(area)’ 


It is in the development of the forcing function that the solution to the momentum 
equation comes into play: the "hat" velocities contribute to the source term of the 
pressure equation. 

In the incompressible case, the second and fourth lines of the above equation disappear 
because the linearization defined in equation (7.85) is unnecessary. The second line is 
treated with the same advection routines that are used for the momentum equation. 

The final step is the velocity update. After the solution for pressure equation, the known 
pressures are used to evaluate the pressure gradients. In order to ensure that a velocity 
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field exists whieh eonserves mass, the pressure term is added baek into the "hat" 
veloeities: 

For SIMPLEF algorithm; 

Equation 7.116. 



Equation 7.117. 


V 


y 




m' 


d(vol)®[Pf 


Equation 7.118. 


Vz 




m' 


cl(vol)'[Pf 


For SIMPLEN algorithm; 


Vx =Vx-- 


— if 

II . .y ^ 
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r« I 'J 


W^|d(vol)"[Pf 


Vy =Vy -■ 


3:1 J 


1 

ay 


f 

i^i J 


V ■ u 

J 


W- 


d(vol)®[Pf 


139 




The global iterative proeedure is summarized below. 

• Formulate and solve equation approximately 

• Formulate and solve “ equation approximately 

• Formulate and solve equation approximately 

• Formulate pressure equation using , and ^ 

• Solve pressure equation for P 

• Update veloeities based on ''^x , ^ and P 

• Formulate and solve energy equation for T 

• Solve speeies transport equations 

• Update temperature dependent properties 

• Solve turbulenee equations for k and ^ 

• Update effeetive properties based on turbulenee solution 

• Cheek rate of ehange of the solution (eonvergenee monitors) 

• End of global iteration 
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APPENDIX D. EXPERIMENTAL DETAILS 


EXPERIMENTAL RIG OVERVIEW 



This photograph shows the test section installed between the plexiglass inlet duct 
and the exit chamber. Note the insulation on the chamber to minimize hat loss to the 
room. The white PVC pipe extending from the exit chamber leads to the flowmeter (not 
pictured here). Although not visible, the thermocouples measuring the outlet air 
temperature are located at the junction of the exit chamber and the PVC pipe. Behind the 
apparatus (from left to right) are the relay board, relay controller and data acquisition 
computer. The inclined manometer can be seen at the lower edge of this photo. 
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TEST SECTION DETAIL 


This photo shows the test seetion diseonneeted from the inlet and outlet duets. To 
the right of the test seetion is the relay board. On this board are mounted the 20 
independent relays to energize the strip heaters. During normal operation this board is 
covered by a plexiglass cover to prevent inadvertent contact with the electrical leads. 
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THERMOCOUPLE PLACEMENT AT TEST ARRAY OUTLET 



This photo shows the placement of the four thermocouples used to find the 
average outlet temperature. The thermocouples are taped to the gasket and placed 
between the exit duct and PVC pipe that leads to the flowmeter. Although the flow is 
expected to be well mixed by this point care was taken to capture various points in the 
flow. As such the thermocouples are placed to measure four different radii (R) in the 
pipe entrance. As seen in the photo, these locations are approximately equal to R, R/2, 
R/4 and R/8. 


FLOWMETER INTERIOR VIEW 



This photo shows the interior of the turbine flowmeter. 
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REYNOLDS NUMBER CONTROL 



These bleed valves were used in conjunction with inlet filters to control the mass 
fiow rate through the test array. In this way the blower (below the lab table) was 
operated at constant velocity for all runs. This provided consistency and repeatability 
between runs. 


144 





APPENDIX E. UNCERTAINTY ANALYSIS 


The uncertainty analysis was conducted using the following general governing 
equation taken from Ref. 21: 



where i? is a given function of the independent variables xi,X 2 , ...x„ 
and wr is the uncertainty. 

Using this method, the uncertainty was calculated for Reynolds number, Nusselt 
number and friction factor. 

REYNOLDS NUMBER UNCERTAINTY CALCULATION: 

Given that Re^,^ = , 

M 


+ ^ 


(note that Wu is assumed to be zero) 
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. ^ r 

Volumetric flow rate: V = — where Vf is the flowmeter voltage output and kf is 

kf 

the flowmeter constant (pulses/ft^). 

From the manufacturer’s calibration sheet, Wv= 0.03 volts and Wk= 1.2 pulses/ft 
andkf = 132 pulses/ft^. Therefore: 



Average area: A = 



where L is the array length 




where x, y and z are dimensions (in mm) in all three axes. 


_ 3 
L 500 


0.006 


_ 'j 'j 

so, (0.017) +(0.006) =0.018 

therefore, 



fMl +(0,018) 

_V J 



For hydraulic diameter: and 
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Where 


= 0.007 
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1/2 
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1 T y 



boOy 
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^Vopen 

V 

open J 


r... \ 


-|l/2 


+ 


W 


Aw 


V y 


(0.017)"+(0.007)' 


1/2 


= 0.018 


Substituting: 


w, 


Ren 


Re 


Dh 




\p) 


+ 


/ \2 
yU j 


+ 




Dh 


y^h j 


1/2 


^0.03^ 


(0.0017) + ^ +(0.018) 




ni/2 


So the relative Reynolds number uneertainty deereases as flowmeter voltage 
inereases eorresponding to higher Reynolds number values. The following table shows 
the range of Reynolds numbers and assoeiated uneertainty values. 


Reoh 

Flowmeter voltage 

Relative uneertainty 

5,000 

0.66 

4.9% 

45,000 

6.30 

1.9% 


NUSSELT NUMBER UNCERTAINTY CALCULATIONS: 


Given Nuj^^ = — 


Nu 


Dh 


^ w'^ 


y n y 


+ 


W 


Dh 


V A y 


1/2 


sinee k is eonsidered eonstant. 


Also, h = ——— and — = 


AAA h 


r V 
vSy 


+ 


w. 


V A y 


+ 


l+A 


V ^ A y 


1/2 


So uneertainty values must be found for Q, A„ and ATim. 

Sinee Q is based on the time that the strip heaters are energized, the relative 
uneertainty is approximately equal to the smallest inerement of “on” time, or 8.547 
seeonds, divided by total run time. For all experimental runs, total run time was 20 
minutes. 
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Wn 8.547 

Therefore, — =-= 0.007 

Q 1200 


For uncertainty in temperature measurements, both the 6 and 12 degree 
temperature difference cases must be considered. With an advertised uncertainty of 0.5K 
for each thermocouple, the relative uncertainty is approximately 0.083 for the 6 degree 
case and 0.042 for the 12 degree case. 


Since these values are much larger than the relative uncertainty values for heat 
flux and wetted area, the relative uncertainty for heat transfer coefficient can be reduced 


^1, 

to: — = 

h 


kQj 


+ 


w 


Aw 


V y 


w. 


V y 


1/2 


nV2 


w. 


V y 


The same argument holds true for Nusselt number uncertainty so, for the 6 and 12 


degree temperature difference cases. 


w,. 


Nu 


-= 8.3 and 4.2 % respectively. 


Dh 


FRICTION FACTOR UNCERTAINTY CALCULATIONS: 

F„r/=^, 

pU^L 


Wf 


f 


w 


A 


Ap 

Ap j 


w 


Dh 


\^h j 




\P J 


\L J 


1/2 


In the above equation, all uncertainties have been found while calculating 
uncertainties for Reynolds number and Nusselt number with the exception of pressure 
drop. Using the inclined manometer, pressure drop uncertainty is constant and is 
approximately 1.25 Pa as the column was graduated in 0.01 inches of water. Therefore, 
the relative uncertainty for pressure drop is dependent on absolute pressure drop as 
shown in the following equation: 


w. 


f 


w 


A 


Ap 

Ap y 


w 


Dh 


J 




\P ) 


\ L y 


n'/2 r 


^\25Pa^^ 


-|l/2 


Ap 


+ (0.018)'+(0.0017)" + 


v500y 
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The following table provides relative uncertainty values for friction factor for 
various pin configurations at each Reynolds number limit. 


Pin type 

ReDh 

Ap (Pa) 

Relative uncertainty 

33 mm round 

5000 - 45000 

40 - 2500 

3.3 - 0.62% 

16.5 mm round 

5000 - 45000 

38 - 840 

3.4 - 0.64% 

10 mm round 

5000 - 45000 

6.4 - 480 

20.0 - 0.68% 

airfoil 

5000 - 45000 

4.5 - 200 

25 - 0.78% 
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APPENDIX F. SAMPLE APPLICATION 


WETTED AREA MAGNIFICATION 

One way of characterizing pin-fin heat exchanger utility is by thinking of the 
array as a wetted area multiplier. Consider a requirement to transfer heat across a 
specific area via forced convection such as in an electronic component cooling 
application. Without a pin-fin array, this would dictate a specific heat transfer coefficient 
that may not be attainable depending on the heat flux requirement. By employing a pin- 
fin array, the wetted area can be magnified to the extent that a reasonable average heat 
transfer coefficient would satisfy the cooling requirement. Geometrically, this wetted area 
magnification factor can be resolved on a unit cell basis. 


First, consider the definition of area density or a: 

^ _ '^welted 

V 

total 

and = 4 ,,, *H 

where is the heat exchanger “footprint” with a specified heat flux 

requirement and H is the heat exchanger height. 


Thus the magnification factor would be; 


A 

aH = 


A 


'face 
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AIRFOIL-SHAPED PIN-FIN CONFIGURATION 

The diagram below shows a unit cell for an airfoil-shaped pin-fin array. Note that 
there are two pins in the cell. It is “5”’ wide, “2A” long and “//” high. 



Using the airfoil dimensions shown below, it is possible to derive an expression 
for the area magnification factor, aH . Although this is derived for an airfoil-shaped pin 
array, the relationship can also be applied to a circular pin-fin array by letting Ry = 0. 
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Wetted area ealeulation; 


A =2* A 

^wetted ^endwall 


- 4 * A A 

^ ^pf ^ ^ ^pinsides 


where 


= 2XS 


pin footprint, = IR^O -IR^R^ 

( R 1 

where 6 = tan ' ^ 

(note that for a eireular pin, yf^^reduees to nR^) 
and the surfaee area of the pin sides, = AHR6 

(for the eireular pin, reduees to ItiR ) 

Thus, 


4XS-4i2*R^0-RR ] + 2(4HR0) 

aH = -i^- - 

2XS 

But, LR = ^ , = ^ = —^,and = 

X D. R-R^, sin 6' 

Substituting, 


-AR(AR^-l) 


aH = 2-- 


sin 0 


sin 6* 




2AR 


LR 


f ^ ^ 


V^>V 


where D^, = D for the airfoil-shaped pin-fin 


and ^ = tan ' 


D^.(AR^-\) 
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While this equation looks imposing, it ean be used to evaluate the trends in 
magnifieation factor when S/D, H/D, LR and AR are varied. Table 17 shows the changes 
in area magnification factor as these basic parameters are changed. Magnification factor 
increases as S/D and LR are decreases and H/D and AR are increased. 

Table 17 Magnification factor trends, airfoil-shaped pins 


S/D 

H/D 

LR 

AR 

ccH 

2 

1 

0.7 

6.2 

2.24 

5 

1 

0.7 

6.2 

2.10 

5 

4 

0.7 

6.2 

2.95 

5 

4 

0.5 

6.2 

2.68 

5 

4 

0.9 

6.2 

3.22 

5 

4 

0.9 

2.4 

3.35 

3 

4 

0.7 

4.24 

3.40 

2 

4 

0.9 

2.4 

5.38 

2 

16 

0.9 

2.4 

13.85 

2 

32 

0.9 

2.4 

26.17 


CIRCULAR PIN-FIN CONFIGURATION 

For the circular pin case, with Ry = 0, 


aH = 2 + 


n 


{X\D){SlD)yD 


H 


2y 


Here, it can be seen that the magnification factor can be increased by increasing 
H/D and decreasing X/D and S/D. However, as shown previously, better performance in 
terms of h vs. E can be achieved with wide spacing pin ratios. Therefore, the best 
practical solution will likely be a tradeoff between desired heat transfer and allowable 
pressure drop. By programming the magnification factor equation and Nusselt number 
and friction factor correlations into a spreadsheet such a trade study could be 
accomplished for a given application. 
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SAMPLE APPLICATION 

Consider an example where the heat flux requirement is 1.6 W/em with a 40 
degree log mean temperature differenee. This defines an effeetive heat transfer 
eoeffieient or heff, requirement of 400 W/m K. 


The required heat transfer eoeffieient, h = 


aH 


Three different pin spaeings and 2 pin height ratios ean be eompared to evaluate 
the best eonfiguration for this applieation. These six eombinations are presented in Table 
18. The fifth eolumn in the table shows the required heat transfer eoeffieient based on 
the effeets of the area magnifieation faetor. The next eolumn shows resulting energy loss 
for eaeh array for a unit faee area as ealeulated from the previously developed 
eorrelations for Nusselt number and frietion faetor. The final eolumn shows heat transfer 
eoeffieient divided by speeifie fluid frietion energy thereby providing a measure of array 
effieieney. 


Table 18 Possible sample applieation solutions 


X/D 

S/D 

H/D 

aH 

h (W/m^K) 

E (W/m^) 

h/E {UK) 

1.5 

1.5 

1 

2.70 

148.3 

29.2 

5.1 

3 

3 

1 

2.17 

183.9 

44.4 

4.1 

5 

5 

1 

2.06 

193.9 

26.6 

7.3 

1.5 

1.5 

4 

6.89 

58.1 

3.1 

19.0 

3 

3 

4 

3.22 

124.2 

37.2 

3.3 

5 

5 

4 

2.44 

163.9 

60.8 

2.7 


Due to the effeets of the magnifieation faetor, the eonfiguration AD =1.5, S/D = 
1.5, H/D = 4.0 would be the best solution for this applieation out of the options 
eonsidered. By way of eomparison, an airfoil-shaped pin-fin array with S/D = 3, H/D = 
4, LR = 0.7 and AR = 4.24 met this eooling requirement with an h/E = 48.8. This 
highlights onee again the superior performanee of the airfoil-shaped pin array. 
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APPENDIX G. LIMITS ON AXIAL PITCH 


SIZE LIMITATIONS BASED ON KNUDSEN NUMBER 

Previously in this work, the minimum heat exehanger hydraulic diameter was 
limited by concerns regarding particulate contamination. To avoid clogging, minimum 
passage dimensions were set at 0.3 mm or 300 microns. This determined the dimensions 
for the optimum practical circular pin case, S/D = X/D = 5.0, H/D = 0.25 andX= 6 mm, 
resulting in a hydraulic diameter equal to 0.6 mm. In a highly filtered environment where 
particle sizes could be on the order of 1 micron or smaller, such as in electronic 
component cooling applications in clean room environments, the heat exchanger could be 
further reduced in scale to improve performance. However, in such a situation flow under 
rarefied conditions of a gas could become a limiting factor. This topic will be discussed 
briefly in this Appendix. The following information on the role of the Knudsen number 
and the Mach number in determining rarefied flow regimes is taken from Ref 25. and 
Ref 26. 

KNUDSEN NUMBER 

The Knudsen number, Kn, is defined as the ratio of the mean free path of the gas 
and the flow characteristic diameter. Thus: 

Kn = — 

A 

A is the mean free path and is defined by: 
kT 

where k is Boltzmann’s constant, 1.380662e-23 J/K 

T = temperature (K) 

p = pressure (N/m ) and 

cr^ = collision diameter of the gas (m) 

For T= 300K,/» = 101350 N/m^ and cr^ taken to be 3 x 10"'*’ m, /I «le-7 m 
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For values ofirn< 0.001, the eontinuum hypothesis is considered valid and the 
flow can be analyzed using Navier-Stokes equations with no-slip boundary conditions. 
For Kn > 0.001, rarefaction effects start to become significant and the flow enters the slip 
regime. In this range the continuum solution validity begins to breakdown. To remain in 
the continuum regime in the current study, the minimum hydraulic diameter should be 
0.1 mm or 100 microns. 


APPLIED LIMIT 

As shown above, the minimum hydraulic diameter must be 0.1 mm to avoid 
rarefaction effects of flow such as slip at the rigid boundaries. In terms of the optimum 
circular pin-fin array, mentioned above, an axial pitch of 1 mm would yield a hydraulic 
diameter of 0.1 mm. This reduction in hydraulic diameter from 6 to 1 mm would result 
in an increase in area density from 65 to 328 cm /cm . 


MACH NUMBER CONSIDERATIONS 

In addition to using Knudsen number as a flow regime check, Mach number can 
be used to further define the region of gas dynamics. From Ref 26, one accepted 


boundary for large Reynolds number flows is determined by the ratio 


Ma 


. If this ratio 


exceeds 0.01, then a slip-flow analysis must be used. Using various inlet velocities, this 
ratio can be evaluated for several pin-fin configurations to see if the slip-flow regime is 
encountered. Table 19 shows various circular pin and airfoil-shaped pin arrays each with 
a hydraulic diameter of 100 microns. In each case, Reynolds number has been 


maximized until 


Ma 


sound is taken as 350 m/s. 


= 0.01. Using an inlet flow temperature of 300K, the speed of 


The Reynolds numbers shown indicate a limit below which the slip-flow region 
would be encountered. Note that especially for the circular pin arrays, the limiting 
Reynolds number is higher than what would normally be considered practical based on 
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the high average flow velocity. From this table it is clear that Mach number 
considerations pose no additional limit beyond the Knudsen number limits described 
above. 


Table 19 Flow regime boundaries based on Mach number 



Circular Pins {X/D = S/D 

= 5.0) 


H/D 

X{m) 

Reoh 

U (m/s) 

Ma 


0.25 

l.OOe-3 

500 

80.9 

0.01 

16 

3.20e-5 

500 

79.2 

0.01 


Airfoil-shaped Pins (S/D = 2.0, LR 

= 0.5, AR = 6.2) 






Ma 

H/D 

X(m) 

Reoh 

U (m/s) 


0.25 

2.70e-3 

45 

7.0 

0.01 

16 

2.30e-4 

45 

7.0 

0.01 
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